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EXECUTIVE SUMMARY 


J.L. Volakis, L.C. Kempel, R.J. Sliva, H. Wang and A. Woo 


Project Goal 

The goal of this project was to develop analysis codes for computing 
the scattering and radiation of antennas on cylindrically and doubly 
conformal platforms. Available techniques and codes can only treat 
antennas on planar surfaces whereas in most cases (missile and 
aircraft platforms), the printed antennas are situated in cylindrical of 
doubly curved surfaces. Modern vehicle platforms are also composite 
and/or coated with dielectric materials and this presents an 
additional challege in the analysis. 

Progress Summary 

Because of our previous experience in finite element and boundary 
integral methods, a hybrid finite element-boundary integral(FE-BI) 
method was employed for the analysis of patch antennas on planar 
and cylindrical platforms. The main challenge in this implementation 
was the efficient evaluation of the cylinder's Green’s function in a 
form useful for the FE-BI method. After 4 months of preparatory 
work, both the FE-BI formulation and the Green’s function evaluation 
were completed and programmed in a rather general purpose 
computer code. This is the first code of its kind for the analysis of 
antennas on cylindrical platforms and because of its future utility, 
several user-oriented features were included with it. The code is 
referred to as FEMA-CYL (Finite Element Method for Antennas on 
Cylinders) and has already been delivered to NASA and Navy users 
along with a test case and users manual. To validate the code, 
measurements were performed by J. Silva at the Naval Weapons 
Center, China Lake and these overlaid the calculations. 

Over the last few months of the project period, we concentrated on 
the extension of the FEMA-CYL code to antennas on coated and 
possibly doubly conformal platforms. Such a requirement precludes 
use of the boundary integral method for terminating the mesh 
because integral equations are inefficient for modeling doubly 
conformal surfaces. Consequently, we examined a class of new 
absorbing boundary conditions(ABCs) for terminating the finite 
element mesh. These ABCs proved quite accurate for antenna 
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applications and their performance and capability is described in the 
first section of this final report. Also, an earlier report describes a 
comparison of the computational advantages associated with ABC and 
boundary mesh terminations(see U-M Radiation Laboratory report 
031 173-1 -T). The implementation of the finite element-ABC 
formulation resulted in the code FEMA-CYLA and we are currently 
preparing a users manual for this code. Note that this code is capable 
of analyzing antennas embedded in coated platforms and results are 
included at least for one of those situations. In addtion, we began the 
development of a finite element formulation for modeling non- 
rectangular printed antennas on doubly conformal platforms. This 
formulation employs prismatic elements which are most suitable for 
coupling with triangular surface meshes in generating the volume 
mesh. For typical patch antennas we anticipate that the mesh 
generation will be done without a need to use a external meshing 
facilities. 

A total of three reports and 4 journal papers were published as a 
result of this consortium agreement: 

REPORTS 

"A comparative study of an ABC and an artificial absorber for 
truncating fintie element meshes" Univ. of Michigan Radiation 
Laboratory Technical Report 03 1 173-1 -T 

This report compares three different techniques for 
terminating the finite element mesh in modeling conformal 
patch antennas and slot arrays. 15pp. 

"Radiation and scattering from cylindrically conformal printed 
antennas" Univ. of Michigan Radiation Laboratory Technical Report 
03 1 1 73-2-T, 166pp. 

This is an extensive report desribing the theory and 
capabilities of the FEMA-CYL code. Many results for scattering 
and radiation by cylindrically conformal antennas are included 
along with mesurements. 

"Radiation and scattering from antennas in coated conformal 
platforms" Univ. of Michigan Radiation Laboratory Technical Report 
031 173-3-T (this report) 


2 



This is the final report and contains several sections which 
highlight the research activity over the length of the 
project/consortium agreement. The report contains the 
following sections: 

1. FE-ABC formulation for patch antennas on a coated circular 
cylinder-- describes the formulation of the FEMA-CYLA code. 

2. Radiation by cavity-backed antennas on a circular cylinder- - 
describes the performance of the FEMA-CYL code for antenna 
radiation analysis and includes the measured results collected 
by J. Silva at the Naval Weapons Center, China Lake, CA. This 
section has been submitted as a paper to /EE Proceedings-Pt. H 
The paper is co-authored with J. Silva 

3. Scattering by cavity-backed antennas on a circular cylinder- 
-describes the performance of the FEMA-CYL code for 
computing antenna scattering. This section is a reprint of the 
journal paper which appeared in IEEE Trans. Antennas & 

P ropagat Vol. 42, Sept. 1994, pp. 1268-1279. 

4. A Hybrid finite element-boundary integral method for the 
analysis of cavity-backed antennas of arbitrary shape- - 
describes the basics of the finite element-boundary integral 
formulation for modeling conformal antennas. This section is a 
reprint of the journal paper which appeared in IEEE Trans. 
Antennas & Propagat. Vol. 42, Sept. 1994, pp. 1233-1242. It is 
co-authored with A. Woo (NASA-Ames) and H. Wang(Naval 
Weapons Center) 


Journal Papers 

J. Gong, J.L. Volakis, A.C. Woo and H.T. G. Wang, "A hybrid finite 
element method for the analysis of cavity-backed antennas of 
arbitrary shape," IEEE Trans Antennas Propagat., vol. 42, Sept. 
1994, pp. 1233-1242. 

L.C. Kempel and J. L. Volakis, "Scattering by cavity-backed 
antennas on a circular cylinder," IEEE Trans. Antennas & 
Propagat., Vol. 42, Sept. 1994, pp. 1268-1279. 
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L.C. Kempel, J.L., Volakis and J. Silva, "Radiation by cavity- 
backed antennas on a circular cylinder," submitted for 
publication IEE Proceedings, Pt. H. 

T. Ozdemir and J.L. Volakis, "A comparative study of an ABC 
and an artificial absorber for truncating finite element meshes," 
Radio Science, to appear in 1995 
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FE-ABC Formulation for Patch Antennas on a Circular 

Cylinder 

Leo C. Kempel and John L. Volakis 
The University of Michigan 
Radiation Laboratory 
1301 Beal Ave. 

Ann Arbor, MI 48109-2122 


Abstract 

The finite element-boundary integral (FE-BI) method has been shown to 
accurately model the scattering and radiation of cavity-backed patch anten- 
nas. Unfortunately, extension of this rigorous technique to coated or doubly 
curved platforms is cumbersome and inefficient. An alternative approximate 
approach is to employ an absorbing boundary condition for terminating the 
finite element mesh thus avoiding use of a Green’s function. In this report, 
a FE-ABC method is used to calculate the radar cross section (RCS) and 
radiation pattern of a cavity-backed patch antenna which is recessed within 
a metallic surface. It is shown that this approach is accurate for RCS and 
antenna pattern calculations with an ABC surface displaced as little as 0.3A 
from the cavity aperture. These patch antennas may have a dielectric overlay 
which may also be modelled with this technique. 


1 Introduction 

Recently, a Finite Element- Boundary Integral (FE-BI) formulation was pro- 
posed by the authors [1] for modeling the scattering and radiation of cavity- 
backed patch antennas recessed in a cylindrical platform. The use of the 
boundary integral for terminating the FE mesh renders the FE-BI method 
numerically exact but leads to a partially full and partially sparse matrix. To 
obtain a fully sparse system, we must use approximate local boundary con- 
ditions for terminating the FE mesh and this is usually done by employing 
absorbing boundary conditions (ABCs). 

In this report, a new conformal ABC recently introduced by Chatterjee 
and Yolakis [2] will be used for scattering and radiation parameter calcu- 
lations in connection with cavity-backed antenna elements on a cylindrical 
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platform. This second order conformal ABC allows the closure surface to 
be brought quite close to the cavity aperture. As a result, the required de- 
grees of freedom are significantly less than the number typically used with 
traditional spherical surface/boundary condition such as the one proposed 
by Peterson [3] or the one introduced by Webb and Kanellopoulos [4], The 
accuracy of this FE-ABC method will be established along with guidelines 
for the distance between the structure and the ABC boundary. In addition, 
this new FE-ABC approach will be used to compute the radiation pattern 
and input impedance of conformal patch antennas with a dielectric overlay. 

2 Formulation 

Consider the computational domain shown in Figure 1. There are two volume 



Figure 1: Typical coated cavity-backed patch antenna with ABC mesh ter- 
mination. 

regions: an exterior region, V 1 , which includes any radome overlay and an 
interior region, V 11 . Both regions may be inhomogeneous and are separated 
by the aperture surface, S ap , and the surface metallization surface, S am both 
of which lie on the surface of the metallic cylinder ( p = a). Thus, the 
exterior region is defined by p > a while the interior region has p < a. 
The computational domain is bounded by the union of the metallic surface, 
S' ' =; S s,n + S C1 " where S cm is the metallic walls of the cavity and the 

ABC surface, S abc . 
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Within the computational volume, the total electric fields may be written 


as 


E{r) = E I {r) + E^'if) reV 1 

= E u (r) reV 11 (1) 

where E cyl (r) = E'(r) + E T (r) as before. The total magnetic fields are 
likewise written 

H(f) = H I (r) + H^if) feV 1 

= H u (r) f € V n (2) 

where H^ij) = H*(r) + H T {r). The boundary conditions are readily written 
in terms of the electric and magnetic fields. Within the cavity, the tangential 
electric field vanishes on the metallic walls 

h x E u (r) = 0 fer (3) 

while on the aperture, the total tangential fields are continuous 

hxE\f) = hxE n {r) r e S ap 

h x H\r) = h x H n (r) - h x f <E S ap (4) 

On metallic surfaces, all tangential electric fields vanish, i.e. 

hxE^if) = hxE^f) = hxE n (r) = 0 r € S sm (5) 

while n x E^if) also vanishes over the aperture 

hxE^if) = 0 fes ap (6) 

since it contains both the incident and reflected fields. Thus, the only non- 
zero electric fields on the surface of the metallic cylinder correspond to the 
unknown fields within each region which are continuous across the surface 
aperture as implied by (4). 

The FE equations may be developed by considering the inhomogeneous 

vector wave equation. Employing the method of weighted residuals and 
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Green’s first vector identity, the weak form of the vector wave equation 
f V x E(p,4>,z) V x Wi(p,<f>,z) 

/ — - — r p dp dq> dz 

Jv, p T {p,<f>,z) 

~ k ° Jv ' Wi{p,4>,z)pdpd<j>dz 

—jk 0 Z 0 l n(p , <t>, z) x H(p, <f>, z) ■ W,(p , <*>, z)dS = f', nt 

J S t 


( 7 ) 


where h(p, (/>, 2 ) indicates the outward pointing normal of the element surface 
associated with the tth unknown, S, is the surface area of that element, and 
H(p, <f>,z) is the total magnetic field. It can be shown that the surface integral 
of (7) vanishes for all elements which do not border the cavity aperture. 
Furthermore, their non-zero contribution is limited to the portion of their 
surface which coincides with the aperture. The interior source functional, 
f- nt , is once again given by 


/; n< = -/ 

Jv t 


J V7 v 

M'(p,(f>,z) 

\ V A 

Pr(P,<t>,z) 


+ J 




z)pdpd(f> <(B) 


where M‘ and J' are impressed currents within the cavity. 

A domain decomposition is accomplished by substituting the total field 
relationships (1) and (2) into (7) and after some manipulation we get the 
FE-ABC equation 

, [ Vxg-Vxft g.g 1 

/ [*» *"■*»*_ 

JV," p r ' 

Is ... (" X V x £') . W, dS = ft 1 + jf., (p x H "") . W, dS - 


Jv u 


1 

pT 


E* 1 ■ Wi dV - jk 0 Z 0 — [ h ■ [Wi X /M dS 

Pre] JS « 1 J 


( 9 ) 


where denotes integration of the aperture associated with the i th un- 
known, S- d is associated with integration over the bounding surface of any 
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dielectric elements in the exterior region, n ri is the relative permeability of 
the elements interior to Sj d whereas /i re is the permeability exterior to that 
surface. 

This set of FE-ABC equations may be written as a linear system of equa- 
tions 



where the FE matrix -4{1— ny may be written as a sum of the FE matrix 
used in the FE-BI formulation [1] and a second term attributed to the ABC 
surface 

^{1-HuJ = M + K U ) (11) 

The new FE-ABC equation (9) is comparable to (7) except that the lat- 
ter utilizes a total field formulation throughout the computational domain. 
However, previously we utilized an integral expression for the total magnetic 
field across the aperture which resulted in the FE-BI equation used in [1] . 
Such an integral expression provides an exact relationship between the total 
tangential electric and magnetic fields over the aperture surface which also 
formed the computational domain boundary. Alternatively, we may employ 
an approximate relationship between these two fields with the goal of re- 
taining the sparsity of the resulting linear system. Additionally, as shown in 
(9), this FE-ABC formulation may be used for coated as well as uncoated 
geometries. In the next section, we will develop an approximate relationship 
suitable for mesh closure. 

3 Conformal ABCs 

Traditional three-dimensional vector ABCs [3, 4] require a spherical outer 
boundary which results in an excessive number of unknowns. New conformal 
ABCs have recently been proposed by Chatterjee and Volakis [2] which have 
an outer boundary that follows the contour of the enclosed geometry resulting 
in a minimal number of unknowns. In this section, the specific expressions 
required by this new ABC for a cylindrical-rectangular box boundary will be 
derived. A definition of ABC order will be given and subsequently the first 
and second-order ABC expressions will be presented. 
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For the purposes of discussion, we define a secondary field as the field 
which is a consequence of equivalent currents that are supported by some pri- 
mary source which is either external or internal to the computational domain. 
Thus for scattering problems, the scattered field is the secondary field while 
the incident and reflected fields are considered primary fields. Likewise, for a 
radiation problem, the radiated field is the secondary field whereas the source 
field due to an impressed current is the primary field. In (9), we recognize 
that an ABC must supply a relationship between the tangential components 
of magnetic and electric secondary fields on the absorbing boundary, S a6c . 

The secondary field may be expressed as a Wilcox expansion 


E a (n,t u t 2 ) 


0 -jk o n 


, ^ lim T 

4 *\fu 


Ei(t u t 2 ) 


up 


( 12 ) 


where u — y/RiRi, Ri — pi + n and pi is a principal radius of curvature. 
In this form, the curvature of the non-spherical wavefront is explicitly used. 
The point of observation is given in Dupin coordinates as 


x = nh x 0 (t\,t 2 ) 


(13) 


where h is the unit normal and x 0 (ti,t 2 ) denotes the surface of the reference 
phase front and therefore, t\ and t 2 denote tangential coordinates on that sur- 
face. Absorbing boundary conditions annihilate outward propagating waves 
up to a certain order. A zeroth-order (P = 0) ABC represents the usual 
Sommerfeld radiation condition. A first-order ABC (P = 1) annihilates all 
fields with up to a u -1 dependency while all higher order fields are reflected 
back into the computational domain. For a cylindrical surface, u = y/p, thus 
the zeroth-order ABC is simply the geometrical optics spread factor while 
the first order ABC annihilates fields up to 0(p - °°). Evidently, as the ABC 
order increases, the reflected fields have an increasingly higher attenuation 
factor and hence the boundary may be placed closer to the geometry without 
inducing erroneous reflections. 

We present the second order conformal ABCs attributed to Chatterjee 
and Volakis[2j. In particular, the appropriate expressions for a cylindrical- 
rectangular box boundary will be given. 

Absorbing boundary conditions provide a local relationship between the 
electric field and its curl which may be approximated as 

nxVxE* = si;+|'Vx [n(n-Vx£ 4 )] +f • V, (n • £ a ) (14) 
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where V< denotes the tangential surface gradient operator. Unfortunately, 
use of (14) would result in an asymmetric system L 4 H Uj due to the last 
term which possesses only one differential operator. An asymmetric system 
requires an iterative solver which utilizes two vector-matrix products such 
as the conjugate gradient squared (CSG) solver such as the one presented 
by [ 5 ]. A symmetric system requires only one matrix- vector product if the 
BiCG solver is used. Additionally, for symmetric systems, only the upper or 
lower triangle of the matrix need be computed and stored. 

The gradient in (14) may be approximated by 

V ( (V-^ S ) = jk 0 V t (h ■ E 3 ) + O (\" v ) (15) 

With both a gradient and a divergence operator present, one operator can 
be transferred to the test vector while the other may remain with the source 
vector. Hence, the resulting matrix may be symmetric since both the test 
and source fields are differentiated. With (15), (14) may be written 

hxVxE 3 = &-E s t +H- V x [n(n- V x £ 5 )] + f - V t (V • £ t s ) (16) 


For the basis vectors used in [ 1 ], V • E 3 is always zero on S abc and hence 
the third term of (16) will not contribute to this form of the ABC. For 
surfaces with a common constant curvature for both tangential directions 
on a surface, this new vector ABC (16) will lead to a symmetric FE system 
However, if the principal curvatures on a surface are unequal, the 
system will be asymmetric. For either (14) or (16), the three coefficient dyads 
axe given by 


a 


0 

1 


s{ z>-a!'-4 

- — - 2/c,l M.) 

J J 


Star — Ak — 2/c,) L 


jk 0 + 3 k t 


(17) 


where K n = Ak = k } — k 7 and D = 4 - 5k^, — In the case of 

(16), 7 must be divided by jk Q due to (15). 
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It is advantageous to consider the second-order ABC for singly curved and 
flat surfaces separately. For a singly curved surface, the unit normal direction 
and curvature parameters are n = p, = — - p and k 2 = 0 as before. After 
some manipulation, we find that (17) becomes 


a 

1 

7 


]2kJTT[7 + D { 3ko + ~ P 

2 p ii , 2 p 

dxp A zz 

j2k 0 p -(- 1 j2k 0 p - 3 


** + , 2 ^- 3 y* +lk ° D 


where D = jk a — Note that /? is not symmetric unless p — * oo 

For the second-order ABC, it is advantageous to segment the matrix entry 
M into three parts. After some vector manipulation, these contributions 
are given by 


j(l)abc 


dS 


j(2)abc __ 

4* i* • v * *i 

h- V x (Jl-W,)] dS 


j(3)abc 


[n ■ W,)\ iS 

(19) 


where /( 3 ) a6c is only used for the original, asymmetric ABC (14). Closed form 
expressions for these three integrals may be readily found by utilizing the vec- 
tor basis functions given in [1]. As mentioned previously, since the principal 
radii of curvature are not identical for a cylindrically curved surface, I^ abc 
will not be symmetric. However, these terms are asymptotically identical as 
the radius of the ABC surface becomes large since the surface will then be 
approximately planar. A symmetric ABC may be obtained by dividing both 
the numerator and the denominator of /( 2 ) afcc by the ABC radius, p, which 
results in a symmetric operator. 


4 Radiation Integral 

Not surprisingly, since the FE-ABC method is being proposed for modelling 
more complex geometries than the FE-BI method presented previously by the 
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authors [1], some support task such as computation of the far-zone radiated 
fields must become necessarily more complex. 

The radiation integral for computing the far-zone radiated fields previ- 
ously was confined to the surface of the cylinder for uncoated antennas. This 
field was solely due to radiating magnetic currents in the aperture of the 
cavity due to the use of a second kind dyadic Green’s function (see [1] for 
details). However, when a overlay or a protruding element is present, the 
radiation integral must contain such material. Thus, for this work, the ra- 
diation integration surface is deformed to contain any material above the 
cavity. To do so, both electric and magnetic currents over this blister must 
be used along with both the first and second kind dyadic Green’s functions. 


5 Results 

The aforementioned FE-ABC has been implemented and in this section, it 
will be used to examine the radiation and input impedance properties of con- 
formal patch antennas with dielectric overlays. However, the formulation’s 
accuracy must first be established via comparison with the FE-BI method 
presented previously by the authors [1]. 

Consider a 2 cm x 3 cm patch antenna residing atop a 5 cm x 6 cm x 
0.07478 cm dielectric filled cavity. This substrate has a dielectric constant 
of 2.17 and the cylinder radius is 15.27887 cm. The patch antenna is fed as 
to excite a pure axial mode (<f> a = 0 °,z a = —0.375 cm). Figure 2 compares 
the FE-ABC and FE-BI formulations. In this, data was taken from 3.0 GHz 
to 3.2 GHz every 5 MHz. The agreement is quite good and we find that the 
resonant frequency is 3.11 GHz. At this frequency, the H- plane radiation 
pattern is shown in figure 3. 

The previous example involved a single antenna element placed within 
a discrete cavity with a minimal ABC surface. Since no superstrate was 
present, the exterior region of the computational domain was limited to the 
immediate vicinity of the cavity. However, if a dielectric coating is used, a 
continuous wraparound exterior region is desirable. In this way, the physics 
involved in substrate mode guided waves is included explicitly in the finite 
element analysis. The next example utilizes such a computational domain 
where the exterior region extends a full 360° around the cylinder while its 
axial lengtii is limited. Furthermore, the discrete cavity which contains the 
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antenna element is minimized. 

To consider the effect of a dielectric overlay, Ke and Wong [6] examined 
the antenna used by Dahele et. al. in [7] which is 3 cm x 4 cm and fed to 
excite a pure axial mode (</>, = 0 °,z, = 1.0 cm). The substrate is 0.0795 cm 
thick with a dielectric constant of 2.32. Figure 4 compares the H-plane pat- 
tern of this antenna for uncoated and with a 0.3975 cm dielectric cover which 
is identical in material parameters with the substrate. In this, each antenna 
was excited at their resonant frequencies, 3.0 GHz and 2.91 GHz, respectively. 
Note that, as one might expect, there is no change between the antenna pat- 
tern associated with the coated and uncoated antennas. The agreement with 
the corresponding E-plane pattern is shown in figure 5. However, a slight 
beam broadening is observed rather than narrowing as predicted by Ke and 
Wong [6], The differences may likely be due to the finite cavity aperture in 
our calculations whereas in [6] the coating and substrate were assumed to 
cover the entire cylinder. Further communication will be conducted with the 
authors of [6] to clear up this area of disagreement. 

In the previous example, the resonance frequency was seen to shift due 
to the presence of the dielectric overlay. Such a shift is illustrated in figure 
6 which shows the input resistance of a 2 cm x 3 cm patch antenna (see 
figure 2) with different overlay thickness. In this example, the substrate and 
superstrate axe identical material (e r = 2.17). 


6 Comments 

In this report, we have presented a new application for the FE-ABC method 
which has been developed at the Radiation Laboratory: radiation analysis 
of conformal antennas mounted on an infinite cylinder. We have sketched 
the formulation and presented some initial validation results. Currently, this 
new FE-ABC formulation has matched data generated by a FE-BI method 
for uncoated patch antennas. We have presented some preliminary data for 
patch antennas with dielectric overlays and are currently seeking appropriate 
measured data. In this report, only axial polarization was considered. Cir- 
cumferential polarization should also be studied with emphasis on dielectric 
overlay effects. 

Initial results are promising and correspond to exportation. It should be 
stressed that this new FE-ABC approach is considerably more flexible than 
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previously developed FE-BI methods since it permits material or protruding 
elements in the exterior region of the cylinder. We will explore the utility 
of this approach for inhomogeneous coatings, complex protruding antennas, 
etc. in future reports. 
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Figure 2: Input impedance for the axially polarized patch antenna which is 
2 cm x 3 cm in a 5 cm x 6 cm x 0.07874 cm cavity. The frequency range is 
3.0 to 3.2 GHz with data taken every 5 MHz. 
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Figure 4: H-plane radiation pattern for a patch antenna which is 4 cm x 3 
cm in a 8 cm x 6 cm x 0.0795 cm cavity with different overlay thickness. 
The dielectric constant of the overlay and the substrate is 2.32 and the feed 
point is (<p a = 0 ^, 2 , = 1.0 cm). 
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Figure 5: E-plane radiation pattern for a patch antenna which is 4 cm x 3 
cm in a 8 cm x 6 cm x 0.0795 cm cavity with different overlay thickness. 
The dielectric constant of the overlay and the substrate is 2.32 and the feed 
point is (< f> s = 0°. z s — 1.0 cm). 
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Figure 6: Input resistance for a patch antenna which is 2 cm x 3 cm in 
a 5 cm x 6 cm x 0.07874 cm cavity with different overlay thickness. The 
dielectric constant of the overlay and the substrate is 2.17 and the feed point 
is (</>, = 0°,2 a = —0.375 cm). 
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Abstract 

Conformal antenna arrays are popular for deployment on curved 
aircraft, spacecraft and land vehicle platforms due to their inherent 
low weight, cost and drag properties. However, to date there has 
been a dearth of rigorous analytical and numerical solutions to aid the 
designer. In fact, it has been common practice to use limited mea- 
surements! and planar approximations in designing such non-planar 
antennas. In this paper, we extend the finite element-boundary inte- 
gral method to radiation by cavity-backed structures recessed in an 
infinite, metallic cylinder. The accuracy of the developed FE-BI code 
for a microstrip patch arrays is established by comparison with mea- 
surements. The formulation is then used to investigate the effect that 
the finite aperture has on the radiation pattern. In addition, the ef- 
fect of curvature on resonant frequency, gain, input impedance and 
pattern shape is examined. 
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1 Introduction 


Modern aircraft and missile designs seek to utilize conformal antenna arrays 
rather than conventional protruding antennas due to their low weight, low 
drag, low cost and flexibility. Although most useful aircraft surfaces possess 
some curvature, the vast majority of available design information assumes 
planar elements. Indeed, the literature is rich with approximate [1], nu- 
merical [2] and experimental [3] design and characterization data for planar 
structures. The most common antenna element is a microstrip patch printed 
on a dielectric coated groundplane. Dielectric coated cylinders have also been 
investigated using approximate [4] and numerical [5] approaches. 

Often, it is desirable to enclose each radiating element within a metallic 
cavity to suppress parasitic substrate coupling [6]. Approximate methods, 
such as the cavity model [1], do not include finite aperture effects since the 
radiating currents are restricted to the immediate vicinity of the patch. Most 
integral equation formulations such as the one proposed by Pozar and Voda 
[2] utilize a grounded slab Green’s function in their construction which pre- 
cludes practical finite aperture simulations. Recently, Aberle [6] proposed an 
integral equation formulation which partitions the geometry into an exterior 
half space and an interior homogeneously filled cavity. This approach deter- 
mines the electric field attributed to the patch and feed currents and the two 
regions are coupled by enforcing field continuity across the finite aperture. 
Unfortunately, as with all integral formulations, the linear system requires 
0{N 2 ) storage and considerable computational effort due the fully populated 
system and the slowly converging cavity Green’s function. The simulation 
of large finite cavity-backed arrays using such an approach is therefore not 
practical. 

An alternative formulation, utilizing the Finite Element- Boundary Inte- 
gral (FE-BI) method, was proposed by Jin and Volakis [7]. This approach is 
also suitable for inhomogeneously filled cavity-backed antennas recessed in a 
metallic ground plane. As with all partial differential equation formulations, 
this approach is associated with a highly sparse system which requires only 
O(N) storage. Additionally, when coupled with a Conjugate or Biconjugate 
Gradient-Fast Fourier Transform (BiCG-FFT) solver, the computational bur- 
den is significantly reduced. The FE-BI method has been successfully used 
for scattering and antenna performance analysis involving planar platforms. 

A similar FE-BI method was proposed by the authors [8] for scattering 
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bv cvlindrical-rectangular and wraparound patch antennas. New divergence 
free, high fidelity edge-based elements were presented along with an efficient 
solution strategy which exploited an asymptotic evaluation of the appropriate 
dyadic Green's function as well as a BiCG-FFT solver. The resulting com- 
puter code was shown to accurately compute the scattering by both planar 
and highly curved elements. This paper investigates the accuracy and utility 
of such a FE-BI formulation for antenna performance analysis. The radiation 
pattern of a single element as well as that of a wraparound array is compared 
with measured results. The importance of modelling finite apertures is exam- 
ined and in addition, the effect that curvature has on the resonant frequency, 
gain, driving point impedance and pattern shape is quantified. 


2 Formulation 


In this section, the FE-BI formulation appropriate for radiation analysis is 
developed for cavity-backed antennas recessed in an infinite metallic cylinder 
(see figure 1). As usual, the finite element formulation permits substantial 
modeling flexibility, including cavity inhomogeneities, lumped loads and mi- 
crostrip feeding lines. 

The FE-BI formulation begins with the weak form of the vector wave 
equation followed by specification of appropriate vector shape functions and 
dyadic Green’s function. The resulting FE-BI equations are then used to 
solve for the total electric fields within the cavity and on the aperture (for 
further details, see Volakis et al [9]). The weak form of the wave equation 
can be written as 




IV, 


VxIV j (p,d),z)-Vxiy,(p,(A,z) 






p dp d(f> dz 


+(k 0 a) 2 6 a (j)S a (i) j s J s [W i (Q,<t>,z)‘p(a,<f>,z)x 

(7 e2 (a, 4>, z) x p(a, <f>, z) ■ Wj(a, 4> , 2 )] d<f> dz d<j>dz j = /‘ nt (1) 


In this, Wi are vector basis functions with support limited to the finite el- 
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Figure 1: Illustration of a typical cavity-backed antenna situated on a metal- 
lic cylinder and the associated coordinate system. 
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ement volume V ■ which is associated with the i th degree of freedom, and in 
a similar fashion. S, and Sj represent aperture surfaces associated with the 
i th and j th degrees of freedom, respectively. The appropriate dyadic Green's 
function is denoted by G e 2 and it has convolutional (© = <? — <p , z = z — z ) 
form when evaluated on the surface of the cylinder, p = a. The unprimed 
coordinates represent the test point while the primed ones denote the source 
point and (p, <p, z) are the usual cylindrical coordinates. The free-space prop- 
agation constant is given by k 0 = jt, where A 0 is the free-space wavelength. 
The cavity is filled with an inhomogeneous material having relative consti- 
tutive properties e r and p r . The function <5 a (i)<$ a (j) identifies when both 
the source and test unknowns belong to the aperture and accordingly con- 
tribute to the boundary integral sub-matrix. The FE-BI equation (1) may 
be rewritten in matrix form as 



[C] [ 0 ] 
[ 0 ] [ 0 ] 


E 7 \ = {0} 

E] nt ) { fl nt } 


( 2 ) 


where the entries of [.4] are due to the FE portion of the formulation and [Q\ 
is the boundary integral sub-matrix. In (2), Ej P and E™ 1 denote degrees of 
freedom associated with the aperture and interior fields, respectively. In (2), 
f\ nt represents the internal excitation and for this paper, a radially oriented 
probe feed is considered. 

The vector elements, dyadic Green’s function evaluation, matrices [^4] 
and \Q ] and the far-zone field formulae are given in [8] and are therefore not 
repeated here. The interior source function is given by 


f,nt = 


-L 


J n Y 

'M in w,T 

1 

Pr{p,d>,z) 


1 


+ jk 0 Z 0 J mt (p,4>,z) } • Wi(p,4>,z)pdpd4>dz 

(3) 


where J ,nt and M' nt are the impressed electric or magnetic current densities. 
For a radially (p) directed probe feed, the impressed current located at (<£,, z a ) 
is given by 

J' nt = p/o—— ————————— (4) 

P 
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and (3) becomes 


fr = 


-A^ | n(g[( 0> -0, 5,)] 


(5) 


if the edge-based elements introduced in [8] are used. 

Having computed the finite element and boundary integral matrices as 
well as the internal excitation, the BiCG method is used to solve for the 
unknown electric fields throughout the computation domain. Given that the 
FE matrix is highly sparse and symmetric, the BiCG method is a well-suited 
choice among iterative solvers. It is also important to note that the matrix- 
vector product associated with the boundary integral can be performed using 
FFTs resulting in a reduced storage and computational burden. The com- 
puted electric field within and on the aperture of the cavity may now be used 
to compute antenna parameters such as the gain and the input impedance. 

The radiated magnetic field is computed by integrating the aperture fields 
with the far-zone dyadic Green’s function given in [8] 

H r (r,6,(t>) = jY 0 k 0 a I G e 2{r,9,(f>-,a,<f) , z') ■ 

J s 

[p(a, <t> , z) x £(a, <f> , 2 ')] d<f> dz (6) 

where (r ,0,0) indicates the observation point in spherical coordinates. The 
far-zone electric field may be obtained from (6) in the usual manner. The 
radiation pattern, directive gain and other useful antenna parameters may 
be calculated using (6). For example, the antenna gain may be computed 
from the far-zone electric field as 


<?«(#,*) = 101og 10 [47r|E'(9,<6)| 2 ] +i01o glo Li- 

L L ZtQrtin 


(7) 


where /£„ is the input resistance which is given below and E T is the radiated 
electric field as r — ► 00 . 

In addition to the antenna gain, designers are concerned with the in- 
put impedance of an antenna for feed line matching purposes. The input 
impedance is comprised of two contributions [10] 


Zin — Zp + Zp> 


(8) 
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where the first term is the probe's self-impedance and the second term is 
the contribution due to the presence of the patch. The probe self-impedance 
accounts for the finite radius of the probe and hence can be omitted when a 
zero-thickness probe is assumed. Accordingly, the driving point impedance 
may be found try calculating the voltage between the patch and the cavity 
base 

Zi n — Zq = -Jj f E(p,<p,z) ■ J' nt (p,4>,z)pdpd<i>dz {[)) 

where the impressed current is given by (4), K refers to the volume of the 
finite element which contains the probe-feed, E(p,<f>,z) is the interior field 
and I Q is the constant current impressed upon the probe. Substituting (4) 
into (9), for the cylindrical edge elements presented in [8], it follows that the 
input impedance if given by 


where 


z,„ = £z;, 

t = l 




i n ( Pfc\ 

I 0 Ctihi y p a J 




( 10 ) 


(ii) 


is the contribution due to one of the four radial edges of the element con- 
taining the probe. 


3 Results 

Two types of antenna elements are considered in this paper and they are 
shown in figure 2 where each patch is ao° x 6 in size with a denoting the 
radius of the cylinder. Although the FE-BI method permits mixed-mode 
feeding, for this paper it is convenient to consider only the two lowest order 
non-hybrid modes. A patch whose radiating side walls are axially oriented 
is termed an axially polarized patch and is fed at <f>, — Circumferentially 
(or azimuthally) polarized patches have radiating walls forming constant z- 
surfaces and are typically fed at z s = The 6 = 90° cut is the E-plane 
for circumferentially polarized patches and the H-plane for axially polarized 
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(a) 



(b) 


Figure 2: Illustration of (a) a circumferentially polarized- patch element; and 
(b) an axially polarized patch element. The radius of the cylinder is denoted 
by a. 
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elements. Observation is considered in the 9 — 90° cut/plane since creeping 
wave effects are a primary interest of this paper. 

Several computed and measured antenna patterns have been published 
for patches printed on a coated cylinder. One such patch, which is 3.5 cm x 
3.5 cm, was used by Sohtell [11] to compare the accuracy of the cavity model 
[4] to a surface current integral equation [5]. The measured data was taken at 
2.615 GHz for a metallic cylinder which was 63.5 cm long and had a radius of 
14.95 cm. The cylinder was coated with a 0.3175 cm uniform dielectric having 
relative permittivity of t T = 2.32. Data was taken for —180° < <fi < 180° in 
the 6 = 90° plane. Figure 3 compares these measured patterns with data 
generated using the FE-BI method for an identical patch placed within a 
360° x 7 cm cavity. This wraparound cavity was chosen to simulate the 
coated physical test body. Note that the H-plane patterns are symmetric 
due to the symmetric placement of the feed, whereas the E-plane patterns 
are not symmetric. The placement of the feed was not specified in [11]; 
however, the agreement for the E-plane pattern shown in figure 3 indicates 
that the position used in the FE-BI model (a 4> a = -1 cm) is reasonable. The 
feed was placed at z 3 — -1 cm for the axially polarized (H-plane) case. 

We next consider patch arrays. Being a rigorous method, the FE-BI 
formulation accounts for mutual coupling and cavity termination effects. The 
H-plane pattern of a four element array was measured to gauge the accuracy 
of the FE-BI approach. Each element is 2 cm x 3 cm and placed within a 
5 cm x 6 cm x 0.07874 cm cavity which is filled with a dielectric having 
e r = 2.17. The cylinder is 91.44 cm long and has a radius of 15.24 cm. 
The cavities are placed symmetrically around the cylinder (e.g. a patch is 
centered at 0°, 90°, 180° and 270°). Only the patch centered at 0° was excited 
while the remaining patches were terminated with a 50fl load. The driven 
patch was axially polarized and the feed was located at z a = —0.375 cm. 
Figure 4 illustrates the excellent agreement between the FE-BI formulation 
and the measured data. 

In a previous paper [8], discrete wraparound cavity arrays were found to 
have a significantly lower radar cross section (RCS) compared to a continuous 
wraparound cavity array. Thus, the size of the cavity had a significant effect 
on the scattering properties of the array. The two antennas presented by 
Sohtell [11] were placed within individual cavities which were 7 cm high and 
approximately 30°, 50°, 90°, 180°, 270° or 360° in angular extent. Figure 
5 illustrates that the azimuthal cavity size has little effect on the radiation 
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Angle (<}>) [deg] 


Figure 3: Comparison of measured and computed data for a circumferentially 
polarized element (E-plane) and an axially polarized element (H-plane). The 
antenna (3.5 cm x 3.5 cm ) was printed on a 14.95 cm cylinder with a 0.3175 
cm coating (e r = 2.32). The probe feed was place at (a (f> a ,z s ) = (-1.0, 0.0) 
for the circumferentially polarized patch and at (a <t> a ,z a ) = (0.0, -1.0) for the 
axially polarized antenna. 
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Figure 4: H-plane pattern for a four element patch array. Each patch is 2 
cm x 3 cm and the patches are placed symmetrically around the cylinder. 
Only the patch centered at 0° is fed while the other patches are terminated 
with 500 loads. 
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pattern for a circumferentially polarized element and a similar comparison for 
the axially polarized patch is shown in figure 6. The back lobe of the antenna 
(near o = 180°) is very small for cavities less than 180° but increases for larger 
cavities. For cavities which lie solely on the forward face of the cylinder, the 
substrate modes apparently diffract off the cavity walls; an effect which has 
little influence on the main lobe of the pattern. However, for wraparound 
cavities and cavities which extend into the back side of the cylinder, the 
substrate modes either shed like creeping waves or diffract strongly behind 
the cylinder thus giving rise to the back lobe. 

Having examined the effect of cavity size on the radiation patterns, we 
will now look at the effect of curvature on the gain of patch antennas on 
cylindrical platforms near resonance. Each of the two antennas used above 
were separately placed within 10.5 cm x 10.5 cm cavities which were em- 
bedded in cylinders with increasing radius. The frequency was allowed to 
vary from 2.4 GHz to 2.7 GHz and the gain (7) was recorded every 5 MHz. 
Figure 7 illustrates that the gain decreases with increasing element curvature 
for a circumferentially polarized patch. Since the input impedance is only 
slightly affected by curvature as shown in figure 8 the resulting decrease in 
gain implies a decreased radiated power which is expected since the effective 
aperture area observed normal to the patch is reduced as the curvature in- 
creases. The axially polarized patch exhibits a greater decrease in gain with 
increasing curvature as shown in figure 9. For this polarization, the input 
impedance is affected by curvature as shown in figure 10. The enhanced 
sensitivity of the axially polarized antenna is due to the combined effects 
of a decreased radiated power and driving point impedance. Both antennas 
exhibit a small decrease in resonant frequency (less than 1.5 percent) with 
increasing curvature as illustrated by the rotation of the curves in figures 8 
and 10. 

The radiation pattern of a circumferentially polarized antenna at res- 
onance exhibits reduced creeping wave interactions with decreasing curva- 
ture due to attenuation as shown in figure 11 when excited at a resonant 
frequency. For circumferentially polarized E-plane observation, the radiat- 
ing surface fields are aligned along the <j >- axis which results in little pattern 
broadening. The radiation pattern of the axially polarized antenna broadens 
as the curvature increases which is illustrated in figure 12. Since the radiat- 
ing aperture fields are aligned along the z-axis for axially polarized H-plane 
observation, the pattern exhibits broadening due to the orientation of the 
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Figure 6: Effect of cavity size on the H-plane radiation pattern of an axially 
polarized patch antenna. 
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Figure 7: Gain (7) of a circumferentially polarized patch antenna near reso- 
nance for various cylinder radii. 






Figure 8: Input impedance of a circumferentially polarized patch antenna for 
various cylinder radii. The frequency range was 2.4 GHz to 2.7 GHz and the 
cavity size was 10.5 cm x 10.5 cm. 
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Figure 9: Gain (7) of an axially polarized patch antenna near resonance for 
various cylinder radii. 
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Figure 10: Input impedance of an axially polarized patch antenna for various 
cylinder radii. The frequency range was 2.4 GHz to 2.7 GHz and the cavity 
size was 10.5 cm x 10.5 cm. 
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Figure 11: Variation of the radiation pattern shape with respect to curvature 
for a circumferentially polarized antenna. 
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Figure 12: Variation of the radiation pattern shape with respect to curvature 
for an axially polarized antenna. 
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surface field with respect to the observation direction. 


4 Conclusions 

In this paper, we presented a rigorous analysis of the radiation by individual 
and arrays of patches placed in a cavity recessed in a cylindrical platform. A 
finite element-boundary integral code was developed and the data generated 
by this code for a typical cylindrical-rectangular patches were found to com- 
pare favorably with measurements. Since the cavity model does not include 
mutual coupling and the usual integral equation formulations requires large 
storage and computational resources, the FE-BI formulation is especially 
attractive for array analysis. 

This FE-BI method was used to study the radiation properties of circum- 
ferentially and axially polarized patch antennas. The finite cavity size was 
found to have little effect on the circumferentially polarized E-plane pattern. 
However, for the H-plane pattern of an axially polarized element, the back 
lobe is significantly larger for cavities which extend from the front side to the 
back side of the cylinder. A wraparound antenna exhibited the largest back 
lobe implying that this lobe is a result of creeping wave shedding rather than 
diffraction off the lateral metallic walls of the cavity. The presence of a back 
lobe must be considered when designing low observable, jam-resistant anten- 
nas or antennas on complex platforms (e.g. an antenna near an obstruction). 
Thus, as was the case for scattering reduction, it is advisable to configure 
the patch antenna in the smallest possible cavity. 

The effect of curvature on resonance, gain, radiation pattern shape and 
input impedance was studied. Both circumferentially and axially polarized 
antennas were considered and was found that the resonant frequency in- 
creased with increasing curvature for both antennas. The gain of both types 
of patch antennas decreased with increasing curvature with the axially po- 
larized antenna exhibiting greater sensitivity attributed to the orientation 
of the radiating surface fields and the decreasing driving point impedance. 
The radiation pattern for axially polarized antennas broadens with increasing 
curvature while the corresponding patterns for circumferentially polarized an- 
tennas does not broaden; however, creeping wave interactions are reduced for 
the latter element with decreasing curvature as expected. We conclude that 
axially polarized antennas exhibit more sensitivity to curvature as compared 
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to their circumferentially polarized counterparts. 
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Abstract Conformal arrays are popular antennas for aircraft, 
spacecraft, and land vehicle platforms due to their inherent 
tow weight, drag, and observables. However, to date there has 
b«n a dearth of rigorous analytical or numerical solutions to 
■id the designer. In fact, it has been common practice to use 
limited measurements and planar approximations in designing 
such nonplanar antennas. In this paper, we extend the finite 
element-boundary integral method to scattering by cavity-backed 
structures in an infinite, metallic cylinder. In particular, we 
discuss the formulation specifics, such as weight functions, dyadic 
Oreen s function, implementation details, and particular difficul- 
ties inherent to cylindrical structures. Special care is taken to 
ensure that the resulting computer program has tow memory 
demand and mini m al computational requirements. Scattering 
results are presented and validated as much as possible. 


I. Introduction 

C ONFORMAL antenna arrays are attractive for aircraft, 
spacecraft, and land vehicle applications since these 
antenna systems have low weight, low drag, flexibility, and 
cost advantages over conventional protruding antennas. The 
majority of previous studies pertaining to nonplanar conformal 
antennas has been conducted experimentally due to a dearth of 
rigorous analysis techniques. Traditional rigorous techniques 
involve an integral equation and are limited in terms of radius 
of curvature and structural complexity. Some approximate 
methods have been considered, but these are restricted in 
accuracy and element shape. 

Recently, the finite element-boundary integral (FE-BI) 
method was successfully employed for scattering analysis 
of large cavity-backed planar arrays [1], The resulting system 
is sparse due to the local nature of the finite element method, 
whereas the boundary integral submatrix is fully populated! 
However, by resorting to an iterative solver such as the 
Biconjugate Gradient (BiCG) method, the boundary integral 
subsystem may be cast in circulant form, allowing use of the 
Fast Fourier Transform (FFT) in performing the matrix-vector 
products. This BiCG-FFT solution scheme ensures 0{N) 
memory demand for the entire FE-BI system and minimizes 
the computational requirements. 

In this paper, the FE-BI formulation is extended to scatter- 
ing by aperture antennas conformal to a cylindrical metallic 
surface. In contrast to the planar aperture array, the imple- 
mentation of the cylindricaJly conformal array requires shell- 
shaped elements rather than bricks, and the required external 
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Green s function must satisfy the boundary conditions on the 
surface of the cylinder. In its exact form, this Green s function 
is an infinite senes that imposes unacceptable computational 
burdens on the method. However, for larue-radius evlinders 
a suitable asymptotic formula is available and herein used for 
an efficient evaluation of the Green s function. In addition 
the resulting BI system ,s again cast in circulant form to 
ensure an 0 ( A ) memory demand and to take advantage 

of the FFT s efficiency when carry ing out the matrix-vector 
product. 

A primary difficulty in studying cavitv-backed antennas 
mounted on curved surfaces is the lack of reference data 
In this paper, scattering calculations based on the FE-BI 
method are compared with data based on different techniques. 
Although such validation is necessarily limited, it provides 
confidence in the formulation's accuracy so that this ap- 
proach may be used in extending the available reference 
data. 

An alternative approach for terminating the FE mesh is to 
use an absorbing boundary condition rather than the exact 
boundary integral used herein. The finite element-absorbing 
boundary condition (FE-ABC) method is associated with 
a higher CPU cost because of its enlarged computational 
domain; however, it is more flexible than the FE-BI method 
presented in this paper since it may include a complex 
radome as well as the cavity-backed antenna elements. Such 
an FE-ABC formulation will be the subject of a future 
paper. 


II FE-BI for Circular Cylinders 

In this section, the FE-BI formulation is developed for 
cavities recessed in an infinite metallic cylinder, having walls 
that coincide with constant p-, <p- or c-surfaces (see Fig. 1). 

As usual, the finite element formulation permits substantial 
modeling flexibility, including cavity inhomogeneities, lumped 
loads, super/substrate antenna configurations, or microstnp 
lines and so on. 

The FE-BI approach possesses both low memory and 
computational demand when implemented with a BiCG- 
FFT solver. Although the system of equations associated 
with the FE formulation is sparse, the boundary integral 
submatrix is fully populated. However, if the aperture mesh 
is a uniform grid, the BiCG-FFT solver may be employed 
for that portion of the system, thus retaining Oi X \ memory 
demand for the entire system. In addition, the solver has 
low computational demand since the sparse matrix-vector 
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■CEMPEL AND VOLAKIS SCATTERING BY CAVITY -BACKED ANTENNAS 



>■ 1 • Illustration of a typical cavity-backed antenna situated on a metallic 
inder and the associated coordinate system. 


xlucts require O(N) operations per iteration and the 
>crete convolutions that utilize FFTs require only 0(N\ogN) 
erations per iteration. 

The FE-BI formulation begins with the weak form of the 
:tor wave equation followed by specification of appropriate 
:tor shape functions and a dyadic Green’s function. The 
ulting FE-BI equations are then used to solve for the total 
ctric fields within the cavity and on the aperture (see, for 
unple, Volakis et al. [2]). For the specific configuration at 
id, the weak form of the wave equation can be written 


f j V x ifyp.Q.z) ■ V x W^p.t.z) 
J\, [ p. T {p.0.z) 


^o f r(P- 0 - j[p. 0 . z) ■ W x [p. 0 .z) jpdpdip dz 

+ (k 0 <i) 2 b a {j )b a (t ) j j ,(a. <t>. z) ■ p(a. <t>. z) x 


G 2 (a. 0 .z) x p(a.O . z ) • W } (a .0 

= ir + ir . 


d<t> dz d<t>dz 


( 1 ) 


this, Vf, are vector basis functions with support over 
volume V;. which is associated with the ith degree of 
dom. and in a similar fashion, 5, and Sj represent aperture 
aces associated with the ith and y th degrees of free- 
i. respectively. The appropriate dyadic Green's function 
enoted by G 2 . and it has convolutional (0 = <p - <p ' . 


; = z - z ) form when evaluated on the surface of the 
cylinder p = a The unpnmed coordinates represent the test 
point, while the pnmed ones denote the source point The 
free-space propagation constant is given by k v - ^ . w here 
A 0 is the free-space wavelength. The cavity is filled with an 
inhomogeneous material having relative constitutive properties 
£ r and p. r . The function £ 0 (»)6 0 (j) is the product of two 
Kronecker delta functions. Hence, it identifies which pairs of 
unknowns belong to the aperture and accordingly contribute to 
the boundary integral submatnx. The right-hand side contains 
an internal source (/ 1 * nf ) and an external source ( f; xt i term 
The former is used only for radiation analysis and is omitted 
for this paper. The latter is used for scattering analysis and is 
discussed later. 

The FE-BI equation (1) may be rewritten in matrix form as 


L41 


{£7} 


\[G] [ 0 ] 


[0] [0]J{£7<} 


{£7} 


ur 1 } 

{ 0 } 


( 2 1 


where the entries of [*4] are due to the FE portion of the 
formulation and [Q] is the boundary integral submatnx. In (2), 
E j P E ) nt denote degrees of freedom associated with the 
aperture and interior fields, respectively. 

An important factor in choosing the finite elements for 
gndding the cavity is the element’s suitability for satisfying 
the mathematical requirements of the formulation as well as 
the physical features of the antenna system. Traditional node- 
based finite elements associate the degrees of freedom with 
the nodal fields and have proven unsatisfactory for three- 
dimensional electromagnetics applications since they do not 
correctly represent the null space of the curl operator, and 
hence spurious modes are generated [3]-[4]. In contrast, edge- 
based elements correctly model the curl operator and therefore 
the electromagnetic fields. In addition, edge-based elements 
avoid explicit specification of the fields at comers where 
edge conditions may require a singularity. Jin and Volakis [6] 
presented edge-based brick elements, which are convenient for 
rectangular- type structures and cavities. For cavities residing 
in a circular cylinder, shell elements are the natural choice. 

Cylindrical shell elements possess both geometrical fidelity 
and simplicity for cylindrical-rectangular cavities. Fig. 2 illus- 
trates a typical shell element, which has eight nodes connected 
by twelve edges: four edges aligned along each of the three 
orthogonal directions of the cylindrical coordinate system. 
Each element is associated with twelve vector shape functions 
given by 


^12 {p.<t>.z) = W p (p. 0 .Z 

'• 0 T . Zf. + ) 

^43(p. </>. Z) = \V p (p. 0 . z 

■Oh z,.-l 

^56(p.<J>. 2) = \Vp{p. 0 . Z 

■■ 0T- Zt ■ - ) 

^87 {P-0.Z) = W p (p. 0. Z 


W’u(p.0. z) = w 0 {p. 0.Z 

:P6--.2f. + ). 

^23 (P-0.Z) = Wdp.O.Z 

Pa- 

M 5b(p. 0 . z) = H 8 (p.u.: 

Pb ■ -Cfc. -). 

^67(p.0,z) = WJp.O.Z 

p a . ■■ Zb- +) 
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Fig 2 Cylindrical shell element. 

W\f,{p.4>.z) = W z {p. <t>. z-.p b ,4> r ,-,+), 

^26 iP-4>.z) = W z (p,<t>.Z\p a ,<j) r ,-,-) 
W 4 s(p,<p.z) = W z (p,<f>.z-,p b ,<p tr ,~), 
W 3 t(p,4>.z) = W t (p.<j>,z;p a , <*>,,-.+) (3) 

where W lk is associated with the edge which is delimited 
by local nodes ( l.k ). as shown in Fig. 2. As seen from (3), 
three fundamental vector weight functions are required for the 
complete representation of the shell element. They are 

»■„(>>. i. i. j) ■= *& l*.- *><* - *> j 

ah p y 

= l-( p _ p)( 2 _ 

H' 2 (p,0.2:p.0.i,5) = ^(p-p)(0-0)z (4) 

where the element parameters <p a , p 6# 0,, 0 r , z 6 , z f ) are shown 
in Fig. 2, £ - p 6 - p a , Q - 0 r _ and h = z t - z b . 
Each local edge is distinguished by p, 0 t z, and s as given 
in (3). The ^-term, which appears in the definition of the p- 
directed weight (4), is essential in satisfying the divergence 
free requirement, i.e., so that V W : = o. 1 Note that as 
the radius of the cylinder becomes large, the curvature of 
these elements decreases, resulting in weight functions that 
are functionally similar to the bricks presented by Jin and 
Volakis [6]. Having specified the vector basis functions, we 
may proceed to develop the matrix entries for the system (2). 

The FE-BI system is composed of two parts: a sparse FE 
matrix and a fully populated BI submatrix, as shown in (2). 

1 n j{(>. o. : ) Will satisfy this requirement only within the volume of the 
element These weighting functions introduce artificial charges on the faces of 
the element and are not divergenceless at element interfaces. This is allowable 
since these elements do not guarantee normal field continuity across the 
element faces. 


The FE matrix entries are represented by 


where constant material properties have been assumed within 
each element. The subscripts \i.j) refer to the row and column 
of the matrix entry and correspond to the test and source edges, 
respectively. The auxiliary functions 

* V x H t(p. 0 % z: p t . 0,. z t , $ t )p dp do dz 
I** 1 ~ ^'•{P-0-z:p } .0 } .z J .s } ) 

•W^t(p, 0. r; p t , 0 t . z x .s l )pdpdodz (6) 

are identically zero unless both test and source edges share 
at least one element in common, resulting in a highly sparse 
system. Physically, such a system is a consequence of the 
locality property inherent in a partial differential equation 
formulation. In (6), the direction of the source and test edges 
are represented by (s.t) € {p,0, 2 }, respectively. Since 
the edges of the mesh are aligned along three orthogonal 
directions, only six unique combinations of (s, t ) are required 
for /* 1J , and only three such combinations for I^ 2 \ Since (6) 
is symmetric with respect to source and test edges, the FE 
matrix will also be symmetric. Evaluations of (6) using (4) 
are presented in the Appendix. 

A lumped impedance post may be included in the formula- 
tion by adding a term to (1) and equivalently to (5); surface or 
subsurface metallization layers may also be modeled. Radially 
oriented lumped loads are approximated in the FE-BI formu- 
lation by a filamentary load located at (0 L , z L ) [2]. Such posts 
have length /, cross-sectional area s and impedance Zi, The 
contribution to [*4] is given by 

A>} = jkoZo - L [ 6J*_-<t>L)6(z-ZL) 

sZ L Jv p 

’ Wi{p, 0 , z)Wj(p, 0 , z)pdpd<j>dz ( 7 ) 

which may be readily evaluated in closed form. In addition, 
infinitesimally thin metallization layers may be represented 
by simply fixing a priori the weight coefficients to zero for 
weights associated with edges which are tangential to the 
metal. This is a consequence of using a total electric field 
formulation. The symmetry and sparsity of the FE system [.4] 
is maintained after the addition of these loads, while the BI 
system [ Q ] remains fully populated and symmetric 
The boundary integral provides an exact boundary condition 
for mesh closure, and its construction relies on a cylindrical 
dyadic Green s function. The entries of the boundan integral 
submatrix are 

= (k Q a) 2 I j W t (a, 0 . z; 0 t . £, . *, * 

Js t JSj 

p(u. 0 . z) x G2<a. <p, 2) x p(a. o.z 

• W s (a, 0 . z ; pj. 0j. Zj . Sj )do dz do <1: ( 8) 
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patta 


Fig. 3. Geodesic paths on a circular cylinder 


where the weight functions are given by (4) and evaluated 
at_the surface p = a. In (8), the dyadic Green’s function 
iG?) satisfies both the radiation condition, and the Neumann 
boundary condition at p = a. This dyadic Green* s function 
may be expressed exactly [8] 


G zz (a. 4> s z) 


(2tt) 
G oz (a. cp. z) 


i_ f * r > 

^ °S ^ H n i2) h) e 




t r 

{2lr) n = ?oc J-~ V k l<n ) Hr! 2 ) (-y) 
°*(a.0.z) 

= — y r - 

(2t) 2 J_ -v 


' Hr! 2 ' {-,) 

( nk; \ 2 H^h)' 

. Hiri-,) 

\ko akj // n ( 2 , (-y)_ 


e At *0-k,z) dkz (9) 


where " = k p a and k f , = y/A 2 — A 2 . However, for large radius 
cylinders (e.g., ka > 3). (9) is computationally prohibitive. In 
these cases, which are of main concern in this paper, it is 
advantageous to employ an asymptotic expression for G? [9]- 
[12]. These employ a creeping wave senes expansion of which 
only the two direct path contributions (see Fig. 3) are retained. 

The formula due to Pathak and Wang [9] 

G :z (a. o. z) ~ — - — -ae _ - , *' >, 

2t r * 


• | i cus’b — qi \ - q,\2 - 3r, 

G^.a.o.z - ' 

■ sin ft ros ft / < 1 - 3(/i 1 - q , i • n .( 


G oc i a. o. i i 


l k -...-jk* 


-q f 


(si n 2 S + q{ 1 ~ q){2 - 3sin : * ■ i r .1 


<?i,ser 2 0(u! Ti - r( J))] 


!()► 


where = 
accurate 


COS 2 S ^ ^ 


s/2fc ? aj 1 ~ r. 

In the definition of J, s is the usual geodesic path 
2 ) and 6 is the direction of the 


and q ~ 


has proven quite 


length (s = \jXa$) 2 w _ tv _ 

geodesic trajectory (ft = tan" 1 [ 2 ^]). Depending on which 
of the two direct paths (shown in Fig. 3) is used. <f> = o or 
$ = 2ir — 0. The soft and hard Fock functions. utci) and 
v(i3), respectively, are characteristic of on-surface creeping 
wave interactions and have been extensively investigated by 
Logan [13]. These functions are also presented in the appendix 
of this paper. Although computation of the Green's function 
(10) is now tractable, evaluation of (8) must be done so that a 
discrete convolutional system is maintained and the singulantv 
of (8) at s = 0 is properly treated. 

Care must be taken in evaluating (8) so that the overall 
storage requirement remains O(N) and the singular integrals 
of (8) are accurately computed. If uniform zoning is used, the 
resulting submatrix ([£]) is block Toeplitz and hence amenable 
to solution using the BiCG-FFT method. For the nonself- 
cell contributions, midpoint integration may be used while a 
regularization procedure must be employed for the self-cell. 
Bird [12] noted that (8) recovers the metallic screen Green s 
function when Q — 0 within the available approximation 
order. This suggests that (8) may be regularized by adding 
and subtracting from (10) the function 

R = if- f | (11) 


2 G 0 (a. <p. e) = 


7 VV 

+ JT 


2ttR 


which is the free-space dyadic Green s function multiplied 
by two. The resulting regularized Green's function (curvature 
contribution) is given by 


</! 


G Z2 (a.0.z) ~ (ros 2 ft + v (l 

(2 - 3cos*0)) [<•( J) - 1] J 

(^-qe~ jk "' sin ft cos #|(1 - 37 G - qi ) 

r(J)-l]l 


G° : (a.o.z) - ; 
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( ; 


' ~~ - J —W jk ' I 1 *nr(t — q ; ] - q, 
•i 2 - .Wrtfij n.-*. - i: 

— q\se.c‘8l u( J > — n J p )’ 1 


( 12 ) 


and since u is no longer singular it ma\ be evaluated numeri- 
cal!) The planar contribution mas be calculated in the manner 
described previously b> Jin and Volakis [5], The FE-BI matrix 
has noss been fully developed and it remains to specify the 
excitation Junction for external sources. 

HI. Plane Wave Excitation 

Plane wave excitation of the geometry is considered in this 
section for scattering analysis. The use of the exact boundary 
condition in (1) allows coupling of an exterior excitation 
field into the cavity. We will describe the form of the source 
function /, xr and discuss its numerical implementation. 

The forcing function, due to exterior sources (/"') is given 
bv 


is obtained from traditional modal analysis These expressions 
may be approximated by retaining onK a few terms of 
the senes if k„a sinff, is small However, as this parameter 
becomes large (e.g.. for large a and 8, - 9(ri i. dpi mas 
be replaced with equivalent asymptouc representations similar 
to those considered earlier. Utilizing Watson s transformation 
and Fock theory [14] m connection with ( 16). we find that 

H‘ : v> ~~ — Y 0 sin ■> sin 8,e }k ~ c ° 561 '- 


t 


,-jkpa sin 0, 


t m< P r 


H? , ~.j2Y 0 c<* 


rrr 


pjk 0 COS 0. c 


k Q a sin 0, 


=jZ 0 k 0 a J \Y,(a.o.z 

x H cyi (a . 0 , z )do dz 


p=i 

- Y 0 sin cos 0, e ]k ° cos e ' z 
2 

. 1 ) p e~ J * QQ sm e ' 

P=1 




) • pia. O . c ) 


“ J 


m 


k Q a sin 0 X 


g n> (m<i > p : 


(13) 

where H,(/j.o.:) is the testing weight for the ith row of the 
matrix and H cyl represents the magnetic field on the cylinder's 
surface in the absence of the cavity. A plane wave 

-jk 0 ik : r j 


E‘ = p 
H‘ = Y, 


p(sin ~ cos#, cos o, - cos -> sin o , ) 


o(sin "y cos 8, sin o, + cos y cos ®, ) 


- e sin y sin 8, 


pjk„\a sin 6, coso, cos 6, j 


(14) 


is assumed to be incident on the cylinder from the direction 
where -» is the polarization angle and e' = (('cosy + 
o'sin -v is the electric field polarization. In these, the difference 
between the observation and incidence angles is denoted by 
o, = 0 - 0 ,. The total surface field is given by the sum of 
the incident and corresponding scattered field from the infinite 
metallic cylinder [14], Specifically, 

H ry, (a.o.;) = H' la.o.r) -r z) 

= oH7' + iH c . yl (15) 


(17) 

in which <f>! = f - (o - 0 ,), <t > 2 = ( 0 - c.) _ 2 , 

m = • 31,(1 complex conjugation is denoted by an 

asterisk. The appropriate far*zone Fock functions 
and / (0) ) are given by Logan [13] and are also presented in 
the appendix of this paper. 

The asymptotic formulas (17) are quite accurate except 
m the geometrical optics region (® 0 ,). In this case, 

Goriainov's [15] expressions 

m yl ~ 

- }osinasin^ie- , * < ’ co5e - z |e~ j *° as ‘ n# ' <t ’ 1 * 

sin 0, cos {<$-$>») \r*t i± i wi* 

+ C *'[G(-mcos(0 - 0,))] \ 


H7‘~-j2Y 0 cosa 


m* 


k 0 a sin 0, 


cos t 


where 

H ry,! 


f _jk, ro> fi : 

id.o.z) = -2Y,,— ; V 


+ e’ k ' a sm c ° s ; [F( - m cos (0 - 0, ) )] * 


COS 


L^ 2 (k 0 a sin 0, ) 




sin - cos 0, 


Av.uMi.fr, H,! 2, (k 0 aYm 8, ) 
. sin ^ 


4* Y o sinacos0 1 e jk ° cos *-~ J e~ jk ° a sin ^«*i 


,jr\( j i-o-C: > 




H7ui.o.:\ = j2Y c ,^—e^ cm(t ' : Y 

nk 0 a 


m ,, . 

J ky^7 l 9 ,{rn 




/n . _ ; 


H,i~ ( k u a 0, f 


( 16 ) 


_ sin e t cos(o-o, ) 


TTl 


G(-mcos (0 -<?>,)) 


~ J kj sin 6i G(1, i-m.c<x(o - 0 ,)) 


(18) 
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have been found to be more accurate and can be used instead 
of (17). The Fock functions (G. G 1 and F) are again defined 
in Logan [13] and given in the appendix. These surface field 
expressions may be used to calculate the entries of the column 
vector {/“*} efficiently via a numerical evaluation of (13). In 
particular, the modal senes (16) is used when fc n asin0, < 10 
and either ( 17) or (18) for k n a sin <9, > 10. as appropriate. With 
the excitation function and the FE-BI matnx now specified, 
the BiCG-FFT method [ 16]— [17] mav be used to determine 
the unknown electnc fields within the cavity . 


IV. Scattering 


Once the cavity aperture and volume electnc fields have 
been determined by solving (2) for an external excitation, the 
radar cross section (RCS) may be calculated The far-zone 
fields may be computed by integrating the aperture fields with 
a suitable Green's function. In this section we present the 
relevant formula for calculating the far-zone fields and hence 
the RCS due to excitation by a plane wave (14). 

To determine the far-zone fields, we begin with the integral 
representation for the scattered magnetic field in terms of the 
aperture fields. We have 


H s (r.8.<t>) = jY 0 k 0 a J G 2 (r. 


8.0: a. 0 . z )• 


Jp(a. o . z ) x E(a. <t> , z 


dd> dz 


(19) 


They must therefore be recast in another form by mplov- 
mg Watson's transformation and Fock theory. a> wj> done 
previously (17). In doing so. we obtain 


-180 ^ fo^: 

4 7T 


Q bo ^ x w I c J*q _ i d t ~ J k a v;:i 




L 


V 0, (m4> p ) -j — - u : -m*, 
k, t a sin 8 


Qtiz ^ fcp sin 8 Cu . H; f ~ ,k a sirN'I; _ 
4rr “ - 


2a7r sin 8 


jk^COsHz Mil N4* ; \f-' 

p=i 


\ __ i 


where the Fock functions are the same as those used with 
(17) due to reciprocity. As was the case for the plane wave 
source, Goriainov’s [15] approximations are more accurate in 
the geometrical optics region (<p ^ o) and similar expressions 
may be obtained for (22). as was found for (17). The far-zone 
scattered field can be computed numerically by using (19) and 
either the series or asymptotic formula as appropriate. Having 
done so, the RCS is calculated from 


o{8.<p) 


lim 47 rr 2 

T — O C 


\E*(r. 8. o)| 
| E'{r. 6. o)\ 


(23) 


with (r, 8 , <t>) indicating the observation point in spherical 
coordinates. When the observation point is very far from the 
cylinder, the dyadic Green’s function in (19) can be replaced 
by its far-zone representation 


G 2 (r. 0. (t>: a. o . ; 




G 0 8<t> 4- G” z 8z 4 G**W 


(20) 


where the unprimed unit vectors are functions of the ob- 
servation position and the primed ones are functions of the 
integration point in (19). The components of this far-zone 
Green’s function 


7 2A; 0 cos 6 ^jk 0 cos 8z 

O ~ — ■ — — — : ; — — — c 


(2n) 2 (k 0 a sin 8) 2 


y i — 

n fTc tfn <2 '(Gasin0) 


(0-0 )) 


3 2 


(27 r)= a 


jk Q cos 8z 


1 


p jn( ~ +(0-0 )) 


n T? x Hn mi (k o as\n0) 

1 

(27t) 2 a sin 8 ' 


jk 0 cos 8 z 


Y. 


1 




H'n\k 0 a sin 6/) 


( 21 ) 


are determined by a mode matching procedure. As one might 
expect, these senes converge rather slowly for large k 0 as\n8. 


Above we presented a FE-B1 formulation suitable for mod- 
eling cavity-backed structures embedded in a circular cylinder. 
Next, we consider a few numerical calculations aimed at 
validating this formulation and giving us an appreciation 
on how the cylinder's curvature influences the scattering 
parameters. 


V. Results 

Having solved for the electric fields induced by an incident 
plane wave, the resulting RCS data must be validated with 
known results. As previously mentioned, available measured 
or computed data is rather scarce, and as a consequence we are 
forced to rely on limiting cases in order to validate this work. 
As the radius of curvature decreases, a cylindrical-rectangular 
cavity will approximate a planar-rectangular cavity. Another 
limiting case involves comparison of an elongated 3-D cavity 
with a corresponding 2-D cavity for normal incidence itf, = 
90°). Finally, we may compare our infinite cylinder results 
with a finite body of revolution (BOR) model tor certain 
polarizations and angles of incidence. We begin with the 
quasiplanar case. 

The first validation effort for scattering b> cavity-backed 
patch antennas relies on the fact that a small patch on a very 
large radius cylinder is quasiplanar and approximates rather 
well an equal sized planar patch. For our test we chose as a 
reference a planar 3.67b cm x 2.75 cm patch residing on a 
7.34 cm x 5.334 cm x 0.1448 cm cavity filled with a dielectric 
having c r = 4. The equivalent patch on a 32.6 cm cylinder is 
6.46° x 2.75 cm residing on a I2.9() c x 5.334 cm x 0.1 44> 
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Angie (0) [deg] 

fig. 4. Compansoo of RCS for a planar parch (3.678 cm x 2 75 emi 
nwirlrn^ on a ,.34 cm x 5.334 cm x 0.1448 cm cav,tv filled w.th fr = 4 

a corres P°ndmg quasi-planar patch on a large radius (32 6 cm) 
cylinder. The operating frequency is 9.2 GHz 



Angle (*) [deg] 


4 ° 51,77 fi „ , MoM rcsul,s and re -B' RCS results ,or 

nrfius of 1A ° CaV ' ty tha ' ‘ S reCfSSed J c >''"dcr with 


a 

a 


cm cavity. At the operating frequency of 9.2 GHz. the cylinder 
has an electrical radius of 10A o . 

Fig. 4 shows the results for the patch on a large radius 
cylinder with corresponding data for the planar cavity-backed 
patch. Clearly, the two RCS patterns are in excellent agree- 
ment, and although Fig. 4 illustrates only monostatic scattering 
in the 4> = 0° plane, additional runs for normally incident 
monostatic scattering and various bistatic situations yield 
similar agreement. 

, ,? 0 K T P sT SOnS may 4,50 made for elongated cavities and 
2-0 MoM results. Long narrow cavities have very little axial 
interaction for principal plane (9 = 90°) excitation and 
therefore results based on this formulation should compare 
well with corresponding 2-D data. It is well known that the 
RCS of a 3-D scattering body of length L » A 0 is related 
to the corresponding 2-D scattering of the same cross section 
via the relation 



*3D = 2l -p 
\ ^0 


Observation Angle (0) [deg] 

c F ^de 6 for a C 3A n P XA 0f a the fin ^ COmpu,ed v,a the ^-Bl method and a BOR 

a radius of 1 A , , ' fi " edwraparound cav, '- v « a cylinder with 

a radius of 1A„ thal is excited by a normally incidem H-polanzed (o = 90°) 


Such a comparison is shown in Fig. 5 for monostatic 
scattering by a 45° x 5A x 0.1A cavity recessed in a cylinder 
with a radius of 1A 0 for both principal polarizations. Once 
again the agreement between the two results is excellent 
iius providing a partial validation of the formulation for 
ughiy curved geometries. We remark that similar agreement 
las been observed for bistatic scattering in the 0 = 90° 
)lane. 

The planar approximation eliminates the effects of curva- 
ure, which is a primary interest in this work, and the 2-D 
ompansons done above are only valid for normal incidence 
nd observation. To consider oblique angles and a highly 
urved structure, we resort to comparisons with a body of 
evolution (BOR) code for wraparound cavities. Since the 
OR code can only model finite structures, we simulate an 
•finite cylinder by coherently subtracting the far-zone fields 


of the finite structure without a cavity from similar data 
which includes the cavity. Such an procedure mimics common 
measurement practices and was found suitable for near normal 
incidence and quite acceptable near grazing incidence in the ' 
case of H-polarization (a = 90°). An example calculation for i 
the former case is given in Fig. 6. where a bistatic scattenne 
pattern is presented in the c> = O' plane due to a plane 
wave incident a. (9, = 90'.o, = 0 =,. Clearly, there ,s good 
agreement between the FE-BI results and data based on the 
BOR formulation. 

The previous comparisons serve to validate the formulation. 

Having done so, it is instructive to examine the effect that 
curvature has on the scattering properties of cavity-backed 
patch antennas. Consider a 2 cm x 3 cm patch residing on a 
cm x 6 cm x 0.07874 cm cavity that is filled with a dielectric I 
having £ r - 2.17. The cylinder has a radius of 15.28 cm Figs 
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tg. 7. RCS frequency response for a 2 cm x 3 cm patch residing in a 5 
m x 6 cm x 0.07874 cm cavity with « r = 2.17 as a function of curvature 
>r E-polanzation (a = 0°). 


and 8 illustrate the behavior of this geometry as a function 
f frequency and curvature. 

Evidently, the resonance behavior of this patch is sensitive 
3 curvature for both principal polarizations. The frequency 
ssponse for E-polarization is more sensitive to curvature since 
le radiating surface field component is parallel to the long side 
•f the patch and cavity. If the patch and cavity were oriented 
o that the long side is in the <t> direction, the response to H- 
olarization would exhibit greater sensitivity. Such an effect 
s important to low observable antenna designers since they 
vant to operate the antenna in the region of lowest RCS. This 
ow return region is a consequence of delicate cancellations 
lue to the physical layout of the aperture. Such cancellations 
re not as complete for highly curved structures as they are 
or planar cavities. 

Conformal antenna designers often use wraparound an- 
ennas to achieve omnidirectional coverage. Two different 
onfigurations are typically used: a continuous cavity where 
he cavity is filled with a single continuous collar of dielectric, 
rnd discrete cavities symmetrically placed around the circum- 
erence of the cylinder. These two configurations are shown 
n Fig. 9. 

Since near resonance the radiation properties of these two 
ypes of antennas are similar, any RCS advantage one might 
assess could govern the appropriate choice of arrays. Fig. 10 
rompares the E-polarized monostatic scattering at 3 GHz in 
he 6 = 90° plane for a wraparound cavity and four discrete 
;avities, where the patches and cavities are identical to those 
jsed in the previous example. The radius of the cylinder 
s 15.28 cm, and the four patches are centered at 0°, 90°, 
[80°, and 270°. Not surprisingly, the wraparound structure 
tas a higher return due to coupling within the substrate, 
■iowever, since in this case the scattered field is due to 
he 2 component of the surface field (^-directed magnetic 
:urrents), both cavities yield large scattered fields in the four 
directional lobes. Fig. 11 is the corresponding comparison for 
4 -polarization. 



Fig. 8. RCS frequency response for a 2 cm x 3 cm patch residing in a ’ 
cm x 6 cm x 0.07874 cm cavity with f r = 2 17 as a function of curvature 
for H-polanzauon (o = 90° ). 


In this case, the scattered field is attributed to the o 
component of the surface fields (^-directed magnetic currents). 
Therefore, substrate modes diffract near the patch, resulting 
in discrete lobes for the discrete array while creeping waves 
shed isotropically for the continuous wraparound cavity. Low 
observable designs will favor discrete cavity arrays over 
continuous cavities since the scattering may be channeled 
in preferred directions and the overall scattering level is 
consistently lower. A final example is shown in Fig. 12, where 
we observe that other than the expected higher scattering from 
the wraparound cavity, the scattering behavior of the two 
arrays is very similar. 

VI. Conclusions 

In this paper, we have presented a finite element-boundary 
(FE-BI) integral technique suitable for electromagnetic scatter- 
ing calculations involving cavities embedded within a circular, 
metallic cylinder. This formulation is analogous to the FE- 
BI approach used by Jin and Volakis [1], [5], [6] and may 
accordingly be used for the analysis of scattering by a large 
array of cavity-backed patch antennas. These cavities need not 
be identical, periodically spaced, or homogeneously filled and 
may in fact may possess lumped impedance loads or surface 
metallization layers. The FE approach employs vector finite 
elements that properly represent the electromagnetic fields and 
possess high geometrical fidelity for cylindrical-rectangular 
cavities. Such elements were presented and are analogous to 
the bricks used for modeling rectangular cavities. In addition, 
we presented an efficient method for evaluating the on-surtace 
and far-zone dyadic Green' s functions. The presented tormu 
lation is amenable to solution using the BiCG-FFT method 
provided uniform zoning is used across the aperture, and a^ 
a consequence this implementation has low computational 
and memory demand. We have presented some validation 
of this work with appropnate limiting cases that proxides 
further archival reference data. In addition, we showed h<*w 
this formulation may be used to influence conformal antenna 
designs. 
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Fig 9 niustrauon of rwo types of arrays: (a) continuous wraparound array; (b) discrete wraparound cavity array 



Fig. 10. Comparison of E-polanzed monostatic RCS at 3 GHz for a 
four-patch array placed on a wraparound collar or in four discrete cavities. 
The cylinder radius is 15.28 cm. The patches and cavities are identical to the 
one used in Fig. 7. The observation plane is B = 90°. 



Angie (♦) [deg] 

Fig. 1 1. Comparison of H-polarized monostatic RCS at 3 GHz for a four 
patch array placed on a wraparound collar or in four discrete cavities The 
cylinder radius is 15.28 cm. The patches and cavities are identical to the one 
used in Fig. 8. The observation plane is 0 = 90°. 


APPENDIX; 

Fock Functions 

The asymptotic form of the dyadic Green’s function with 
observation both on the surface of the cylinder and in the 
far field involves Fock functions. These have been extensively 
studied and tabulated by Logan [13]. The numerical evaluation 
of these functions is performed for either small arguments or 
large aiguments. 

The on-surface Fock functions used in this paper are 
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w 2 (t) 


e j(T dr (A- 1 ) 


where w 2 (t) and its derivative w 2 (r) denote Airy functions 
of the Second Kind. For small arguments (£ < 0.6). the 
asymptotic expansion of (A- 1 ) is given by 


~ + 4 ?+ & 
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*t- ... 
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while a rapidly converging residue series is used for £ > 0.G: 
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TABLE U 


j Table B-2 j! 

1 Constants for (A-5) and (A-6) j; 

m 

c(m) 

o (m ) 

.4i(m) 

1 

0.7473831 

1.01879297 

0.5356566 

2 

-0.6862081 

3.2481975 

-0.41901548 

3 

-2.9495325 

4.82009921 

0.38040647 

4 

-3.4827075 

6.16330736 

-0.35790794 

5 

8.9378967 

7.37217726 

0.34230124 

6 

56.1946214 

8.48848673 

-0.33047623 

n 

l 


9.53544905 

0.32102229 

8 


10.52766040 

-0.31318539 

9 


11.47505663 

0.30651729 

10 


12.38478837 

-0.30073083 



12 Comparison of H-polanzed monosumc scattering at 3 GHz by a 
-patch arra> placed on a wraparound collar or in four discrete cavities, 
cylinder radius is 15.28 cm. The patches and cavities are identical to the 
used in Fig. 8 The observation plane is o = 0°. 


TABLE I 


Table B-l 

Zeros of the u’ 2 (t) and w' 2 (t) 
= iT„|e _; 5 and T n = |r'|e- J ? 

n 

w 

l T nl 

1 

2.33811 

1.011879 

2 

4.08795 

3.24819 

3 

5.52056 

4.82010 

4 

6.78661 

6.16331 

5 

7.94413 

7.37218 

6 

9.02265 

8.48849 

7 

10.0402 

9.53545 

8 

11.0085 

10.5277 

9 

11.9300 

11.4751 

10 

12.8288 

12.3848 


10 

u($) ~= 2 e J T v /?^ £ (r„) _1 e“ j€T " (A-3) 

n = 1 


calculated using 

g {0) (O = 2 . 0 e _j4 r £<-1.3 

dm 


= 1.39937+ V — 1.3 < £ < 0.5 

rrv 


10 


kq (m)£ 


= V 0.5 < £ < 4.0 

m= 1 

-(0. 8823-jO. 5094)C-j^ 

= 1.8325e L J £ > 4.0 (A-5) 


s“>({) = -,2.o(e + e < -2-8 

^ y - f ^ Kr 1 - 2.8 < { < o ,5 

" ml 

m= 1 

JiL J*a(m)<] 

05<f£4 '° 

m= 1 

= -1.8325(0.8823 - jO.5094 + ]£ 2 ) 

(0.8823— j0. 5094 )£ -j £ j 


• ei 


£ > 4.0 (A-6) 


ire r n and r n are zeros of wz(t) and u? 2 (t), respectively. 
)se zeros are given in Table I. 

Tie far-zone Fock functions are given by 



f (,) (0 = 
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e^ T 

u ' i ( t ) 


dr 


dr 


G (,) (0 = 9 (l) (Oe ]i r 
F {,) (£) = f in (Oe ]L r 


(A-4) 


;re and its derivative irjfr) denote Airy functions of 

First Kind and the integration contour is given by Logan 
1 These functions, # (0) (£K and / (0) (£), may be 


= 0.77582 -f e~ }w/3 V 

f rn' 


m — l 


C < -1.1 

-!.!<£< 0.5 


10 

- E 

m = l 


e [*a(m)£] 

,4t' (m) 


= 0.0 £ > 4.0 


0.5 < f < 4.0 


(A-7) 


with constant k = e -jW6 and the coefficients for (A-5) and 
(A-6) given in Table II. 

The corresponding constants for (A-7) are given as shown 
in Table III. 
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table m 
Table B-3 


LEL 

I c(m) 

q(tti) 

— 

At (mj J 

| 1 

1 .146730417 

2.33810741 

0.70121082 p 

1 2 

0.86284558 

4.08794944 

•0.80311137 

3 

-2.0192636 

5.52055983 

0.86520403 

4 

-9.977776 

6.78670809 

-0.91085074 1 

5 

-14.59904 

7.94413359 

0.94733571 j 

6 ! 

49.0751 

9.02265085 

-0.97792281 1 

* 

8 

9 

1° 


10.04017434 

11.00852430 

11.93601556 

12.82877675 

1.00437012 
-1.02773869 
1.04872065 
-1.06779386 1 


J 2b 

1) Ss$t 


f )(z — Zt )dz 

J z b 

Mf f Pb 

/ \p- i>s)(p- p t )dp 


J 2 _ *,&th [1 4 1 

lfa)2 4 ( ^ fc Po) * o^» ~Pf\pt - P j I 


2 P,P,{PI -p 2 a ) ,X to - 0,110- Of 1,/c- 


Eacrh of the above unevaJuated integrals is of the form 


FE Matrix Entries 

The matrix entries for the FE portion of the system (5) 
are given in this appendix assuming that the cylindrical shell 
elements (4) are used in (6). These integrals are given by 

/( 1 ) _ [211 ( Pb\ [° r , 

pp ~M7r 6/l n (^J J 0I 

/ (*-■=.)(«- i* )dz 

J z b 

r(i) _ f° r , 

pz q 2 / (0 ~ )dtf> 

uj _ «,5to r /1 

00 (tA)2 [^ 4 ^* -p “) 

+ 3 (^. 4 - p,)(p3 Pi>) + Ip 4 p f ( p 2 -pM 


2(pfc P Q ) 2i(^ a -f p t ) + p,/5, In T — 

\ Pa 


1 .) _ r ( i 1 2 ox / „ 

~ ( ta )- a \ v 2 ~~ + Pt) + PsPf In | — 

V Pa 

1 /*©*• *i 

^ 2 ~ ^ “ 0 a )( & — (t>t )d 0 

''Of 

!' = jjftgt \ n (f>*>\ f°', 

(c*hy 2 l {^ a ) j ci i 0 -^(®-<Pt)d<i> 

I (z - z,){z - z t )dz 
J z 6 

> _ *,5, a [1 1 

(thy- [4 ( pb ~ Pa ^ + 3 f/5s + P'HpI - pI) 

" \p,p,(pI - pi ) j X J \ z - - 2s)(z _ if)d . 




TTie integrals /" M2 ’ are used in the assembly of the FE portion 
<[-4 ) of the system. 

5) 

II 
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A Hybrid Finite Element-Boundary Integral 
Method for the Analysis of Cavity-Backed 
Antennas of Arbitrary Shape 

Jian Gong. Student Member IEEE. John L. Volakis, Senior Member. IEEE . A. C. Woo. and H. T G Wang 


Abstract — An edge-based hybrid finite element-boundary inte- 
gral (FE-BI) formulation using tetrahedral elements is described 
for scattering and radiation analysis of arbitrarily shaped cavity- 
backed patch antennas. By virtue of the finite element method 
(FEMi. the cavity irregularities, the dielectric super/substrate 
inhomogeneities, and the diverse excitation schemes inside the 
cavity may be readily modeled when tetrahedral elements are 
used to discretize the cavity . On the aperture, the volume mesh 
reduces to a triangular grid allowing the modeling of non- 
rectangular patches. Without special handling of the boundary 
integral sy stem, this formulation is typically applicable to cavity - 
backed antenna systems with moderate aperture size. To retain 
an < u V i memory requirement, storage of the full matrix due 
to the boundary integral equation is avoided by resorting to a 
structured triangular aperture grid and taking advantage of the 
integral's convolutional property. If necessary, this is achieved 
by overlaying a structured triangular grid on the unstructured 
triangular grid and relating the edge field coefficients between the 
two grids via two narrow banded transformation matrices. The 
combined linear system of equations is solved via the biconjugate 
gradient (BiCG) method, and the FFT algorithm is incorporated 
to compute the matrix-vector product efficiently, with minimal 
storage requirements. 

I Introduction 

M ICROSTRIP antennas have been extensively investi- 
gated experimentally, analytically, and numerically for 
decades. By and large, numerical methods have been serving 
the engineers and researchers in the analysis and design of 
these conformal antennas for many years. Among them the 
moment method in conjunction with various integral equation 
(IE) formulations played a major role [IJ-[3]. However, IE 
methods are associated with field representations in which the 
appropriate Green's function for the specific geometry must 
be employed, and this limits their versatility. Moreover. IE 
techniques are usually formulated on the assumption of an 
infinite substrate, a model that obviously deviates from the 
practical configuration, leading to inaccuracies for larger band- 
width antennas. Furthermore, in the context of IE methods, 
antenna excitations are represented using simplified models 
that differ more or less from the actual configurations. Also. 
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patch ground plana 


~ Bl aperture area 

c » 

FEM cavity volume ♦ — excitation cable 

Fig 1 Illustration of a cavity-backed patch antenna and the different com- 
putational regions 

due to the singularity of the current distribution near the patch- 
probe junction(s). special measures must be taken [4], not to 
mention additional IE complexities due to possible substrate 
anisotropies or inhomogeneities in the antenna substructure. In 
contrast, the hybrid finite element-boundary integral (FE-B1) 
technique alleviates these difficulties, and this was demon- 
strated recently w'hen the method was applied to rectangular 
patch antennas [5]. 

In this paper, we present an edge-based hybrid FE-BI 
formulation using tetrahedral elements for a characterization 
of arbitrarily shaped cavity-backed antennas. An example of 
such a configuration is shown in Fig. 1. where a cavity is 
recessed in a metallic ground plane enclosing the FEM volume 
and the antenna elements on the aperture may be excited by 
different schemes, such as a simple probe, a magnetic frill 
generator, a practical coaxial cable, rmcrostnp lines, slots, or 
a CPW line In the context of the FEM, the cavity is first 
discretized into a number of tetrahedral elements that naturally 
reduce to mangles on the cavity's aperture. For nonrectangular 
patches this triangular gndding is, in general, nonuniform, 
and the exact boundary integral formulation based upon this 
mesh applies to any patch shape. As a result, the hybrid 
FE-BI technique is capable of modeling arbitrarily shaped 
cavity-backed antenna configurations, different substrate m 
homogeneities, anisotropies, and various practical excitation 
schemes. 

As is w'ell known, the boundary integral (Bli equation 
subsystem leads to a fully populated matrix whose m/c i> 
determined by the number of aperture mesh edges For large 
apertures, this analysis becomes impractical in terms ot Morugc 
and computation time requirements, and to overcome 
inefficiency a uniform zoning of the aperture is requ.rcd 
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By resorting to the structured mesh, the boundary integral 
matrix can be cast into a discrete convolutional form, thus 
permuting the computation of the matrix-vector products via 
the discrete Fourier transform (DFT). avoiding a need to 
store the full BI matrix This memory-saving scheme has 
already been applied to IE solutions involving rectangular 
[6]. [13] and triangular [7] surface grids, and in this paper 
we describe how the BiCG-FFT solver is implemented for 
triangular meshes. The differences between the rectangular and 
triangular meshes are also described, and results are presented 
that demonstrate the method's versatility' in computing the 
scattering and input impedance of various nonrectangular 
printed antennas. 


II. Formulation 

In this section, we present the edge-based hybrid FE-BI 
formulation using variational principles, where the matrix 
algebra notation is employed so that one can readily extend 
the formulae to the general anisotropic case. As derived in 
[5], the complete functional pertinent to the scattering and 
radiation by a cavity-backed configuration (shown in Fig. 1) 
may be written as 

F{E) =\ fJJv{ (V * E ) x E) dv 

+ 2jkoZo J J (E x ET) ■ z dS + E 
• (jk 0 Z 0 J, + V x dv - 2/cq JJ (E x i) 

-{JL tExi) { 1 + k' rv ) G ° iry)dS '} ds 

(i) 

where J, and M, represent interior electric and magnetic 
current sources within the cavity V; H' is the incident field, 
if any. from the exterior region; the surface 5 emcompases 
the cavity aperture excluding the portion occupied by the 
antenna elements; e r and n r denote, respectively, the relative 
permittivity and permeability; k$ is the free space wave 
number, 7 the unit dyad, and Go(r. r') the free space Green’s 
function with r and r' denoting the observation and integration 
points. 

A. Cavity Volume Modeling 

In proceeding with the discretization of (1), it is convenient 
to reexpress it as 


\ i . yl f”' 


with 


= iiTHLHT}:, 


IT \ = 


\F\, = 


(V:\ 


\w 


\eJ 


U — S. V- ; 


in which l Ul is the u ( u = j*. y or r ) component of the volume 
vector basis functions along the z th edge. The unknoun vector 
{ ^ } e h&s six entries, one for each tetrahedron edge iln this 
paper, we use square brackets for matrices and curly brackets 
for vectors.) Inserting (3) into (2), and taking the hrst variation 
of F\ with respect to {£} p , yields 


e 


where 


w -IJL ]f |j^02 0 ^J xy l 


+V x 



[DV}J = 


r ly {K} -£^ 


dr 


<61 


<7 ) 


( 8 ) 


To carry out the above integrations, it remains to introduce the 
volume expansion or shape functions V e . For our implemen- 
tation we employed the linear edge-based shape functions for 
tetrahedral elements given in [8] and [9], 


F = F\ + F s (2) 

where F\ denotes the volume integral contributions and, 
similarly, F s accounts for the surface integral contributions. 
The cavity volume is subdivided into X tetrahedral elements 
C !e = 1.2. ■■■X). and within each tetrahedron the field is 
expanded using edge-based elements as 

£ = ra r {£} f o) 


B. Aperture Modeling 

To discretize the surface integrals in (1). the aperture is 
subdivided into triangular elements since these correspond to 
the faces of the tetrahedrals. Within each triangle, the field is 
represented as 

E = [S]J{E a } r (9, 
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where 





in which 5 U1 is the u(u = r.y) component of the surface 
vector basis functions along the ah edge. On substituting (9) 
into the surface integrals of (2) and taking the first variation 
of F s with respect to we obtain 

= £ {[£],{£.}< + {L},} (11) 


where 



■G 0 (r.r')dSdS' (12) 


and 


{L}r = j2k 0 Z 0 jj s [5 e ] j dS . 


(13) 


Note that in (12) the elements of the array [S e ] are functions 
of the observation vector r, whereas the elements of [S c ] r are 
with respect to the integration point r' . A suitable set of linear 
edge-based surface basis functions is 


S,(r) 


2X ix(r — r -> e(r > 

0 


r € S e 

otherwise. 


(14) 


In this expression. /, denotes the length of the zth edge and r, 
is the position vector of the vertex opposite to the ilh edge. 
Since each edge shares two mangles, one is defined as the plus 
and the other as the minus triangle. Therefore, e(r) is given by 


Hr) 


fi res; 
1-1 res- 


(15) 


w here S, = . The constant A e in ( 14) denotes the area 

of the plus or minus triangle depending on whether r £ Sf 
or r £ Sr ■ We note that S,(r) x £ yields the basis functions 
used bv Rao et al. [10] in their moment method solution of 
boundary integral equations. 


C System Assembly 

To construct the final system for the solution of the electric 
field components we combine (5) and ( 1 1 >. and after assemblv 
we obtain the system 

{\A]{E) + {K}} + {{B}{E 9 } + {L}} =0. (16) 

In this, { K } and {L} are the excitation vectors due to the 
interior current sources and the exterior excitation, respec- 
tively. The unknown electnc field vector {E\ consists of all 
field expansion coefficients with respect to the element edges 
except those coinciding with perfectly electrically conducting 
(PEC) walls, PEC antenna element(s). or PEC pins inside 
the cavity. Finally, the vector {£,} represents the unknown 
surface fields whose entries are part of those in {£} with 
their corresponding edges on the aperture. The explicit ex- 
pressions for the matrices and vectors in (16) can be readily 
extracted from (6), (7), and (12) (see also [1 1]). It is evident 
that [A] and [B] are symmetric as a result of the assumed 
isotropic medium and reciprocity. In addition, [A] exhibits 
high sparsity due to the FEM formulation, whereas [Bj is fully 
populated. Two approaches may be followed in earning out 
the solution of the combined subsystems when an iterative 
solver is employed, such as the biconjugate gradient (BiCG) 
method [13], These two approaches differ in the manner used 
for the evaluation of matrix-vector products called for in the 
iteration steps. One could sum the coefficient matrices [.4] and 
[B] by adding up the corresponding matrix entries pnor to 
the execution of the BiCG algorithm, or instead the resulting 
vectors may be summed after carrying out the individual 
matrix-vector products. We observed that the first approach 
is more efficient in terms of computation time after reordering 
the combined matrix and storing only the nonzero elements. 
This is because, in the context of this scheme, the combination 
of the two matrices is performed only once outside the 
iteration. However, the second approach is compatible with 
the BiCG-FFT scheme, where the FFT algorithm is employed 
to exploit the convolutional property of the integral operator, 
thus eliminating a need to explicitly store the entire BI matrix. 
Below, we discuss the implementation of the matrix-vector 
product of the boundary integral system for the BiCG-FFT 
solution. 

D. Implementation of the Boundary Integral 
Matrix Vector Products Using FFT 

We refer to Fig. 2, which shows an overlay of a uniform 
grid over a nonuniform mesh. The boundary integral equation 
is implemented using the structured triangular gnd, and the 
relation between the unstructured and structured mesh is 
described in the next section. We recognize that the triangular 
grid consists of equal right triangles and thus involves three 
different classes of edges (classes 1, 2, and 3). These include 
the x directed, y directed, and the diagonal edges, all of which 
are uniformly spaced. For the FFT implementation, each class 
of edges is independently numbered in accordance w ith their 
geometric location. Specifically, the /th class will cany the 
numbering (m.n) if the edge is the mth along the x direction 
and the nth along the y direction. The indices (m. n } take the 
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m 0 J 2 3 

F.g 2. Overlay of a structured triangular aperture mesh over an unstrucmred mesh, shown here to conform to a circular patch 


values 


matrix vector product can now be calculated as 


m = 0, 1.2. ■ .A/ 1 * * 

n = 0. 1. 2. • ■ • , A rl 

with i = 1 for the y directed edges, i = 2 for the diagonal 
edges, and ? = 3 for the x directed edges. Consequently, we 
find that 

(TV-1 t = l 
A n = < TV - 1 i = 2 (17) 

[TV -2 i = 3 

where M and TV denote the numbers of elements along the x 
and y directions, respectively. 

To perform the integrations for the evaluation of the bound- 
ary integral matrix elements, it is now convenient to rewrite 
the basis functions (14) in terms of the new indices (m.n). 
We readily find that the edge-based basis functions associated 
with each of the aforementioned class of edges can be rewritten 
as shown in (18M20) at the bottom of the page, where the 
superscripts refer to the edge class. Each entry of the boundary 


[M- 2 i=l 
XV = )\i -i i = 2 
[ M - 1 i = 3 


{B1 subsystem} = [5]{E, } 


-E 


£ £ 


m ' =0 n ^ li 



(21 > 


in which (m.n) are the geometric location indices for the nh 
class observation edges, whereas (m\ n' j are the same for the 
jlh class edges belonging to integration elements. Thus, the 
specification of the indices m. and u completelv defines 
the entry = n\f l 4 A/ of the column resulting after the 
execution of the boundary matrix-vector product. It is readily 
found that 

B' J 

mn ,m n 

— ^°JJs J Is' S™ Tl Sm'n' Go(r.r') dxdyds f dy 

+ (A/A.v)2 jj s jj g . Mr)f ' (r,/ ' / ; 

Go(r. r) dx dy di f dy' (22) 


1 f(nA»-y)i-+(T-TnAx)i ( x.,,}t = *.7 

-rnnt^-vV) - — s {y - ( n 4 l)At/)x 4 ((m 4 2 )Ax - x)y ix.y) e S~ 

otherwise. 


= : lAy}2 


Ax Ay 


J')y € S t 

other 


3 1 j ~ f — * y 1 v ‘* w 

f- 7 * - .V) = ^ (y — n A;</)x 4 (mAx — x)y 


(y - (n 4 1 )A.y)x 4 ( mAx 
0 

((n 4 2)A y - y)x 4 (x - (rn 4 1 )Ax),v f x. y) € S; 
~ r> x ■'* {./■.</] € S; 

otherwise. 


((81 


(19) 


( 20 ) 
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I 


with 


( 1 = 1 

l, = •( v '(Ax) 2 ~(Ay)* « = 2 (23 > 

( Ax i=3. 

More important, it can be shown that the BI subsystem exhibits 
the convolutional property B tJ , , = B' 3 and 

thus we can rewnte (21) as 

3 

[Bl{£,} = Y, 3 ' 1 * E> (24) 

J-l 


where the * denotes convolution. It is now seen that the 
computation of the boundary matrix-vector product can be 
performed by employing the 2-D discrete Fourier transform 
(DFT), thus avoiding a need to store the BI matrix other than 
diose entries that are unique. When the symmetry property of 
is invoked, implying 


B' J , M = B? , 

(m—m , n — n ) (m — ’m.rt — 'Vi) 


(25) 


t is concluded that the total nonredundant entries in the BI 
natrix are 


3 3 

N p = Y, ]T N '( M ' + MJ - !)• 


1=1 1 


(26) 


rhis should be compared to the \ entries whose 

torage would normally be required if the BI system was not 
ast in convolutional form. We remark that N p is nevertheless 
;qual to twice the number of entries required for uniform 
ectangular grids [6] for one class of edges. To avoid aliasing, 
[ is necessary that n _ n ,) = B' p 3 (m,n) be cast in a 

-D array that has the usual periodic form, and zero padding 
nay also be required to make use of the standard FFT routines, 
pecifically. the matrix-vector product (21) is executed by 
sing the MFTxNFT array (27), shown at the bottom of the 
age, with the corresponding field vector given by 


E J p (m.n) 


f E J (m.n ) 0 < m < M J < 0 <n<N J 

\o otherwise. 

(28) 


nd MFT and NFT must be powers of 2 if a radix 2 FFT 
Igorithm is used. 


In the BiCG-FFT algonthm the BI subsystem vector in 
computed as 


3 

{BI subsystem} = ^ 5 { DFT - 1 { DFT { B\; } DFT(£;}} ? 

1 = 1 

(29) 

The presence of the operator S indicates the necessary reorder- 
ing of the 2-D array that results after the inverse FFT operation 
into a single column with the proper indexing for addition to 
the FEM subsystem. It should be remarked that in contrast 
to [7] the integrals (22) are evaluated without introducing 
any approximation. This is necessary' to preserve the global 
combined system symmetry. 

III. Mesh Overlay Scheme for 
Nonrectangular Patches 

As described above, the BiCG-FFT solver requires uniform 
aperture gridding so that the BI subsystem can be put in 
block circulant form. This can always be achieved during mesh 
generation whenever the patches are rectangular in shape or 
in case of radiators that are placed at some distance (usually 
small) below the aperture. However, for circular, triangular, 
or other nonrectangular patches on the aperture, it is not 
possible to construct a uniform mesh using the mesh generator. 
Typically, the aperture mesh is necessary to conform to the 
patch shape, leading to an unstructured free surface grid. In 
this case, to make use of the efficient, low-memory BiCG- 
FFT algorithm, an approacn is proposed to overlay on the 
unstructured aperture grid another coincident structured gnd. 
as shown in Fig. 2. The boundary integral subsystem is then 
constructed using the overlaid uniform gnd, whose edge fields 
can be related to those on the unstructured grid via two sparse 
transformation matrices. That is, it is necessary to append to 
the system (16) the relations 

{£.},. = [I>]{£.} bu 

{£.}»« = {T B ){E a } u ( 30 ) 

where the subscripts u and nu refer to the field coefficients 
of the uniform and nonuniform aperture grids, respectively. 
Also, [TV] and [ Tq ] refer to the forward and backward trans- 
formation matrices, respectively, with N u and denoting 
the numbers of the uniform and nonuniform mesh edges on 
the cavity aperture. 


Bp 7 (m.n) 


B X] (-m. — n). 
B lJ ( -m. -n). 


0 < m < M l 
0 < n < N J 

MFT - M l + 1 < m < MFT 

0 < n < A rj 


) B'Um.n - 1 - NFT), 


0 < m < M % 

NFT - N J + 1 < n < NFT 


B tJ (rh - 1 - MFT, fi — 1 - NFT). 


MFT - A/' + 1 < rn < MFT 
NFT - N 3 + 1 < n < NFT 


(27) 


10 otherwise. 
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Fig 3 Illustration of the parameters and geometry used in constructing the 
transformation matrix elements between the structured and unstructured 

/o T ° d . env ! 11,6 elements of l T r], we begin with the expansion 
(9) and enforce it at three points on each edge belonging to 
the uniform grid. We conveniently place these three points at 
the center and ends of the edge (see Fig. 3). Given the fields 
at these points, we can interpolate the field along the (m.n) 
edge of the uniform grid using the weighted average 


and the global indices , and j correspond to the rth uni; 
grid I edge and the jth nonuniform gnd edge The subscript 
is the global index used in numbering the nonuniform o n _ 
edges, whereas the subscript i ■ = 1. 2. or .1, ,s the local ed-e 
index used ,n the definition of the bas.s functions 5. . se, 
(14)). We remark that the explicit computation and storage 
of the transformation matnx elements results in a substantia 
increase in efficiency because it aso.ds the usual assembh 
process dunng each mteranon step and that the proposed 
overlay scheme allows the analysis of large nonrec.ancular 
patch arrays because storage of a fulls populated BI sv'stem 
matnx is avoided. The user needs only to provide an additional 
data file that flags the uniform gnd edges Iving on a PEC 
element, and this is an important user-onented feature of the 
formulation. Follow, ng the same procedure, we can ohtam 
the expression for the entnes of the backward transformation 
matnx. It should be noted that assummg each uniform and 
e ge traverses three or less nonuniform gnd mangles the 
nonzero entnes in each row of [7>] will be 9 or less However 
they can reach a maximum of 18 if the midpoint and endpoints 
reside on an edge of the nonuniform gnd. The maximum 
nonzero entnes ,n each row of [T B ] will be 15. but the tvp.cal 
number will be much less. 


\ 2A; ndl 

1 

^ rmd 

^mid 

E 

k = 1 

1 

-V. 

_ A 


J «no 2 

+ 2NZZ E ™( r ' nd 2 


ena2 fc=l J 

in which e u denotes the unit vector along x, y or the diagonal, 
epending on the class of edge being considered. The quanti- 
ties nu re P resent the fields in the nonuniform grid triangles 
with the superscript k being a sum variable in case r end 
or r mid specify a point shared by more than one tnangie. 
viously, A endl . A r m i d , and denote the number of 

nonuniform gnd tnangles shanng the node at r endl , r mid . and 
r e nd 2 . respectively, and will typically be equal to unity 
After assembling (31) into (30). we find that the elements 
of the forward transformation matrix are given by 

irr) '' = r- -{wb-i: 

V * = i f=i 

j V m,d 3 

V j y " y " (r m jd ) 

* m,d 

2 3 'I 

^ 2 A ends ^ ( r end 2 ) / (32) 


^ = 1 (=1 


(,JI ~ {o 


1 J = j( 

0 otherwise 


IV. Numerical Considerations 

Based on the presented FE-BI formulation, a computer 
program was written for the analysis of the radiation and 
scattering by cavity-backed patch antennas of arbitrarv shape. 
The antenna geometry is supplied to this program in an 
input file that, as a minimum, must contain lists of 1 , the 
nodes and their ( x.y.z ) coordinates. 2) the nodes formme 
each tetrahedron. 3) the nodes on the cavitv aperture, and 
4) the nodes on metallic boundaries. For arbitrary antenna 
geometries, it is necessary to employ a sophisticated volume 
mesh generation package, and a number of these are available 
commercially. Typically, each of these packages generates a 
universal file” that can be readily preprocessed to extract 
the aforementioned input lists. 

A major effort was devoted to writing the program in a 
manner that minimizes the storage and computational require- 
ments. Specifically, the boundary conditions on the metallic 
surfaces are enforced a priori to obtain a system that involves 
only nonzero field components. The sparse finite element 
matrix was stored as a single array of length A A ,, where 
A, is the total number of unknowns within the cavitv volume 
and An Z denotes the maximum number of nonzero row entnes 
The BI matrix was stored in different ways, depending on 
whether the FFT was to be employed for the evaluation of 
the matnx-vector products. If the B.CG solution , (( be 
earned out without the FFT. then the A. , A. BI imearal 
matrix was added to the FE array, resulting in a l-D amav 
about A,.A nz A; long. For slot antennas, including cavitv- 
backed spirals, and moderately sized svstems. n found 
preferable not to use the FFT. thus avoiding anv , mediation 
errors. In that case the generation of a single ...rrmncj FE- 
BI matnx before execution of the BiCG alg. reduces 
the computational requirements This is he’.uw -umber 
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of operations associated with the repeated combinations of 
the FE and BI subsystems within the BiCG iteration is 
avoided. 

When the FFT is to be used as part of the BiCG solver, 
the FE and BI matrices must be kept in separate arrays 
throughout the execution process. In this case, the FE matrix is 
again stored as a single array, and similarly the nonrcdundant 
elements of the BI matrix are stored in another single array 
of length 9A (2A / -I). The factor of nine is due to the 
three classes of edges, and as usual M and N denote the 
number of elements along the x and y directions, respectively. 
Because of the storage and computational efficiency of the 
BiCG-FFT algorithm, it is necessary to resort to uniform 
aperture grids for conformal antennas involving a substantial 
number of aperture edges. Of course, one should always use 
uniform triangular grids when the patches on the aperture are 
rectangular or if the array supports a superstrate. In the case 
of nonrectangular patches it will be necessary to overlay a 
structured triangular grid over the unstructured grid generated 
by the mesh generator. This must be done in the preprocessing 
stage, and should be taken into account when constructing 
he FE matrix. For scattering computations, the overlay of 
he structured grid is almost always the preferred approach 
because it does not, generally, compromise the accuracy of 
he computed scattering cross section. However, for antenna 
>arameter computations, the interpolation scheme between the 
tructured and unstructured grid edges may be of concern, 
lepending on the specific antenna geometry. Generally, thin 
nnular slots and piannar spirals should be treated without 
esorting to structured gnds, and to our experience this does 
ot cause a large computational burden because these antennas 
re associated with small apertures. In the case of circular, 
iangular, or other nonrectangular patches, the structured grid 
/as not seen to compromise the computational accuracy. Of 
ourse, conclusions based on one type of antenna do not 
ecessarily apply to others, and thus the suggested alternatives 
lust be examined separately for each antenna before choosing 
ne approach over the other. Of importance here is that the 
)mnulation is suitable for modeling any antenna shape and 
;ed structure. 



frequency (GHz) 


Fig. 4. Comparison of the computed and measured back scatter RCS 
as a function of frequency for the shown circular patch The incidence ancle 
was 30° off the ground plane 



V. Results 

We present below some representative numerical results for 
e purpose of validating and demonstrating the robustness 
r the tetrahedral formulation for scattering and radiation by 
fferent configurations of cavity-backed antennas. In each 
se the computed results via the FE— BI method are compared 
ith reference measured or calculated data. 

Scattering and radiation by a circular patch: Fig. 4 illus- 
ites a circular patch residing on the surface of a 0.406 cm 
ick substrate having a relative dielectric constant of e r = 2.9. 
le patch's diameter is 2.6 cm, and the substrate is enclosed 
a circular cavity 6.292 cm wide. This cavity and the patch 
?re recessed in a low cross section body for measuring its 
-S. A comparison of the measured and calculated backscatter 
$ RCS as a function of frequency is also shown in Fig. 4. 
r this computation, the direction of the incident plane wave 


Fig. 5. Comparison of the computed and measured inpul impedance for the 
circular path shown in Fig. 4. The feed was placed 0.8 cm from the center of 
the patch, and the frequency was swept from 3 GHz to 3.8 GHz 

was 60° from normal, and as seen the agreement between 
measurements and calculations is very good throughout 
the 4-9 GHz band. Input impedance measurements and 
calculations for the same patch are displayed in Fig. 5. 
The probe feed in this case was placed 0 cm from the 
patch s center, and it is again seen that the measurements and 
calculations are in good agreement. 

Radiation by a one-arm conical spiral: We considered the 
modeling of this radiator to demonstrate the geometrical 
versatility of the FE-BI method. Two projections of the spiral 
radiator and surface mesh are illustrated in Fig. 6. The top 
and bottom edges of the stnp forming the spiral follow the 
,,TM “ S p - 0.0503Aexp[0.221(<r>±2.66)]. c = «_ expt U.221o 
V>l. PAGE 1 % 
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F.g 6 Illustration of the configuration and mesh of the one-arm conical spiral used for the computation of Fig 7 


where (p.o.z) denote the standard cylindrical coordinates. 
°- ^ e q ual to 0.0832A and 0.0257A. respectively, and 
0 < <p < 2x. This spiral arm resides on an inverted cone (9.24 
cm tall) whose bottom cross section has a diameter of 1.68 
cm and the top cross section has a diameter of 21.78 cm. For 
our calculations A = 30 cm (/ = 1 GHz) and the spiral was 
situated in a circular cavity 10.01 cm deep. The computed E c 
principal plane radiation pattern taken in the <t> = 90° plane, 
using a probe feed at the cavity base, is given in Fig. 7. It 
is seen that this pattern is in good agreement with the data 
given in [12]. As can be expected, the E e pattern (not shown) 
differed from the measured data near the horizon because of 
interference from the finite circular cavity housing the spiral 
that was included in the analytical model. The latter was not 
part of the measurement configuration, which consisted of the 
spiral antenna on a large circular plate. 

Anninar slot impedance: Fig. 8 shows a narrow circular 
(O.i 5 cm wide) annular slot situated in a circular cavity 24 7 
cm wide and 3 cm deep. Because the annular slot is narrow, 
the implementation of the BI subsystem is very small for 
this app lication, and as a result there is no need to invoke 
the FFT in the BiCG algorithm. The FE-BI method is basi- 
cally quite effective in modeling small aperture configurations 
without a need for special computational considerations. Input 
impedance calculations as a function of frequency for this 
radiator, excited by a probe placed across the slot, are shown 
in Fig. 8. and agree well with the values calculated via a 
modal-boundary integral method [14], For these calculations, 
the frequency was swept from 700-1000 MHz. The dielectric 
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6 

Fig. 7. Comparison of the calculated radiation pattern < E j. taken in the 

° = 9 ° C P lane - Wlth data in reference f 1 2] for the one-arm conical spiral 
shown in Fig 6. 

constant of the material filling the cavity was set to * r = i.:r, 
as in [14], and this is an effective value to account for 
the presence of a dielectric slot cover used as part of the 
measurement model for holding the plate. 

VI. Conclusion 

We presented a hybrid finite element-boundary integral 
(FE-BI) formulation that incorporates linear tetrahedrals. The 
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i • 12.35 cm 
b * 0.75 cm 
p 0 ■ 7.7 cm 
0.7 < f < 1 GHz 


d * 3 cm 


ti 


t htt m 


IX) 



Fig. 8. Comparison of input impedance calculations for the illustrated cav- 
ity-backed slot. 


method was specifically developed for the radiation and scat- 
tering analysis of cavity-backed printed antennas, where the 
FEM is used for modeling the cavity region and the BI 
equation acts as a global boundary condition for terminating 
the mesh on the cavity aperture. The FE-BI formulation is 
particularly suited for the analysis of complex configurations, 
and much emphasis was given here in developing a solution 
technique requiring 0 (\) storage in spite of the resulting 
fully populated BI subsystem. The latter was achieved by 
making use of the convolutional property resulting from the 
structured mesh, thus permitting use of the FFT in the BiCG 
iterative solver for computing the matrix-vector products. For 
nonrectangular patch geometry, a novel numerical scheme is 
proposed to overlay on the unstructured free mesh a uniform 
triangular grid, avoiding the storage of the large BI subsys- 
tem. For scattering calculations associated with large aperture 
structures, use of the FFT proved essential in minimizing the 
computational requirements. 

A number of patches, slots, and planar and nonplanar spiral 
antennas were analyzed for the purpose of demonstrating the 
versatility and accuracy of the FE-BI technique. Certainly, 
the need to use a sophisticated mesh generation package is 
deterrent to the application of the technique for the analysis 
of simple antenna configurations. However, this is unavoid- 
able when dealing with complex geometries and. moreover. 


the pervasive use of such commercial packages on desktop 
computers makes the technique quite attractive 
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